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Abstract: In the context of the Minimal Supersymmetric extension of the Standard 
Model (MSSM), we reanalyze the search for the heavier CP-even H and CP-odd A neutral 
Higgs bosons at the LHC in their production in the gluon-fusion mechanism and their 
decays into gauge and lighter h bosons and into top quark pairs. We show that only 
when considering these processes, that one can fully cover the entire parameter space of 
the Higgs sector of the model. Indeed, they are sensitive to the low tan /3 and high Higgs 
mass ranges, complementing the traditional searches for high mass resonances decaying 
into r-lepton pairs which are instead sensitive to the large and moderate tan/3 regions. 
The complementarity of the various channels in the probing of the complete [tan /3, Ma] 
MSSM parameter space at the previous and upcoming phases of the LHC is illustrated 
in a recently proposed simple and model independent approach for the Higgs sector, the 
/iMSSM, that we also refine in this paper. 
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1. Introduction 

The probing of the electroweak symmetry breaking mechanism and the search for possible 
extensions of the Standard Model (SM) of particle physics has become the main mission 
of the CERN Large Hadron Collider (LHC). Among these extensions, Supersymmetry 
(SUSY) [1] is considered as the most appealing one as it addresses several shortcomings of 
the SM, including the problem of the large hierarchy between the Planck and electroweak 
scales. While the search for SUSY was unsuccessful at the first LHC run, the increase of 
the center of mass energy of the machine from 8 TeV to the 14 TeV level will significantly 
improve the sensitivity to the new particles that are predicted by the weak scale theory. 
These consists not only of the superpartners of the known fermions and gauge bosons but, 
also, of the additional Higgs bosons beyond the state with a mass of 125 GeV that has 
been observed by the ATLAS and CMS collaborations in the first LHC phase [2], 

As a matter of fact, in low-energy SUSY scenarios, at least two Higgs doublet fields 
H u and are required to break the electroweak symmetry and to generate the isospin-up 
and down type fermion and the W/Z boson masses. In the simplest scenario, the Minimal 
Supersymmetric Standard Model (MSSM), the spectrum consists of five states [3-5]: two 
charged 17 , a CP-odd A and two CP-even Higgs particles h and H, with h being the state 
observed at the LHC while H is heavier as present LHC data is strongly indicating [6]. 
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The phenomenology of the Higgs sector is described entirely by two input parameters, 
one Higgs mass that is usually taken to be that of the pseudoscalar A boson Ma and the 
ratio tan /3 of the vacuum expectation values of the two doublet fields, which is generally 
assumed to lie in the range 1 < tan (3 < mt/mb ~ 60. This is the case at tree-level where, 
for instance, the lightest h boson mass is an output and is predicted to be M ^ <Mz \ cos 2/3|, 
i.e. Mh < Mz at high tan/3 for which |cos2/3|~l. However, this relation is violated since 
important radiative corrections, that introduce a dependence on many SUSY parameters, 
occur in this context [7-9]. It has been recently shown that, to a good approximation, the 
MSSM Higgs sector can be again parametrised using the two basic inputs tan/3 and M a , 
provided that the crucial LHC information Mh ~ 125 GeV is used [10-12]. 

It is known that two efficient channels can be used to directly search for the heavier 
MSSM Higgs particles at the LHC and probe part of the [tan/3, Ma\ parameter space 1 . The 
first one is the search for light charged Higgs bosons that would emerge from the decays 
of the copiously produced top quarks and would decay almost exclusively into a r lepton 
and its associated neutrino for tan/3 > 1. For almost all values of tan/3, the latest ATLAS 
[14] and CMS [15] results now practically rule out the mass range Mjj± < 160 GeV, which 
approximately corresponds to Ma < 140 GeV in the MSSM. The second efficient channel 
is the search for high mass resonances decaying into r-lepton pairs, which would be the 
signature of the production of the heavy neutral H/A states and their decay into r leptons. 
The rates for this channel can be very large at high tan/3 values, as a consequence of the 
strong enhancement of the H/A couplings to bottom quarks and r-leptons. This process is 
particularly favored as, for a heavy enough A boson, one has the mass degeneracy relation 
Mff ~ Ma that in practice leads to search for a single resonance and allows to combine 
the rates for A and H production. The most recent ATLAS [16] and CMS [17] results 
with the data collected at the first LHC phase, exclude at the 95% confidence level (CL) a 
significant portion of the [tan /3, Ma] plane for sufficiently high tan /3 values. 

Except in the narrow mass range Mz < Ma < 140 GeV, where the lower value corre¬ 
sponds to the exclusion limit from negative Higgs searches at the LEP collider [18,19] and 
the upper one is due to the present limit from charged Higgs boson searches at the LHC 
(which can be straightforwardly interpreted in the [tan /3,M H ±\ parameter space as the 
H± properties depend only on these two parameters in the low mass range), the low tan/3 
region of the MSSM has not been considered so far by the experimental collaborations. The 
reason is that in the benchmark scenarios that are used to interpret the various experimen¬ 
tal limits on the cross sections times branching ratios in the context of the MSSM [20, 21], 
the SUSY-breaking scale is usually set to relatively low values, Ms ~ 1 TeV, that do not 
allow for a heavy enough h state at too low tan /3. Indeed, the radiative corrections to the 
mass Mh depend logarithmically on the scale Ms and, for instance, one cannot obtain a 
value Mh ~ 125 GeV for tan /3 < 3-5 in the MSSM, even if one favorably tunes the other 

1 Of course, there are also indirect limits on the MSSM parameter space, in particular from the measure¬ 
ment of the couplings of the observed h particle at the LHC; see for instance Refs. [10,13]. These limits are 
nevertheless slightly model dependent as, for instance, they can be affected by SUSY particle contributions 
to the h production and decay rates. These indirect limits, that exclude low values of the A, H, H ± masses, 
are complementary to those from the direct Higgs searches on which we will focus in this paper. 
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SUSY parameters that enter the loop radiative corrections, in particular the stop mixing 
parameter Xy which also plays an important role in this context. This is the case of the 
so-called maximal mixing or M™ ax scenario which is defined such that the value of My 
is maximized, i.e. for a stop mixing parameter Xt — \/6Mg in the dimensional reduction 
scheme [20]. The situation is even worse for different values of the Xt parameter. 

In fact, in most of the [tan/3, M 4 ] parameter space, the measured value My, ~ 125 GeV, 
which should be now considered as an essential information on the model, is not satisfied 
in the M“ ax benchmark scenario with Mg = 1 TeV nor in the alternative benchmark 
scenarios that are presently used to interpret the experimental searches in the context 
of the MSSM. If one allows for an uncertainty of a say 3 GeV in the determination of 
the h mass in the MSSM, from unknown higher order contributions for instance [22], the 
situation is acceptable if the h mass is confined in the range 122 GeV < Mj j < 128 GeV. 
Nevertheless, it remains annoying that for each point of the [tan/3, Ma] parameter space, 
one has a different My value in these benchmark scenarios. 

A straightforward and easily implementable solution to this problem has been proposed 
in Refs. [10-12]: if the experimental constraint My ~ 125 GeV is enforced, one in fact 
removes the dependence of the Higgs sector on the dominant radiative correction and, 
hence, on the additional SUSY parameters, in particular Mg and Xt.. One can again 
parametrise the MSSM Higgs sector using only the two basic inputs tan/3 and M/ 4 , exactly 
like it was the case at tree-level. The masses of the heavier H and H± states as well as 
the mixing angle a in the CP-even sector are given by very simple expressions in terms of 
tan/3 and Ma with the constraint My = 125 GeV. It was shown that this approximation is 
very good in most of the MSSM parameter space that is currently accessible at the LHC, 
even when subleading radiative corrections are also considered [ 10 ]. 

In this minimal and almost model independent approach, called the hMSSM in Ref. [10], 
one has access to the entire [tan /3, AL 4 ] parameter space without being in conflict with the 
LHC data, as the information My = 125 GeV is taken into account from the very begin¬ 
ning (this is not always the case for the Higgs couplings which conflict with the measured 
ones at low Ma)- In particular, the low tan/3 region can naturally be accessed, but at the 
expense of assuming a very high SUSY scale Mg. The reason is that at tan/3 values too 
close to unity, the tree-level h mass becomes very small, MyPdMz\ cos2/3|—>-0. To increase 
My to ~ 125 GeV, the radiative corrections that grow logarithmically with Mg need to be 
maximized and hence, a very large scale, Mg > 0(100) TeV for tan/3 < 2, is required. 

The low tan/3 region can be directly probed by the search for the heavier H/A (and 
eventually H±) states and for relatively low Higgs masses, M yy ~ At a <350 GeV, two ways 
have been suggested. First, one can use the same constraint discussed above from the search 
of resonances decaying into r-lepton pairs [12], Indeed, the rates for A/H production are 
appreciable at low tan /3 as the dominant process, the gluon-fusion mechanism, is now 
primarily mediated by loops of top quarks that have significant couplings to the H/A 
bosons; at the same time, the decay of at least the A boson into rr pairs has a still 
appreciable rate. The second way is to reinterpret the existing ATLAS and CMS exclusion 
limits from the search for a heavy SM-like Higgs boson decaying into a pair of massive 
gauge bosons [23, 24] in the context of the MSSM. At low tan /3 and not too large Myy 
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values for which we are not yet in the decoupling regime with a vanishing H coupling to 
massive gauge bosons, the rates for the decays H —> VV with V = W, Z, as well as for 
gg-^H production, are still significant. In addition, searches for the resonant hh [25,26] 
and hZ [27] topologies have been performed at the LHC with the available 25 fb -1 data 
at \fs = 7+8 TeV, and one can reinterpret them in the context of the MSSM where the 
production cross section for gg—^H/A and the branching ratios for the decay modes H-^hh 
and A—^Zh below the tt threshold can be substantial; see Refs. [12,28]. 

The two types of searches mentioned above, with results that were preliminary and ob¬ 
tained with a subset of the LHC data collected at yfs = 7+8 TeV, have been used in Ref. [12] 
to set constraints on the [tan /?, Ma] plane; excluded regions have been delineated using 
some approximations and extrapolations. In the present paper, we update this discussion 
first by using the latest ATLAS and CMS results, especially the final H/A^r t + t~ and 
t —> bH + —> brv analyses [14-17] as well as heavy SM Higgs searches in the H —> WW, ZZ 
channels [23, 24], with the full set of 25 fb -1 data collected in the first LHC phase. In ad¬ 
dition, constraints from more appropriate analyses in the A^-hZ and H^-hh topologies 
where the resonant case has now been considered [25-27] will be included. We will then 
extrapolate these results to estimate the sensitivity of the 14 TeV LHC run, with at least 
an order of magnitude higher integrated luminosity than the one accumulated so far. 

Above the tt threshold, i.e. for Ma,h f 350 GeV, the previously discussed search 
channels will have little relevance at low tan/3 values as, because their couplings to b quarks 
and r leptons are not enhanced anymore, the heavier H and A bosons will dominantly decay 
into tt pairs, the top-quark Yukawa coupling oc mt/ tan (3 becoming then large. As already 
mentioned, the main Higgs production channel will be the gluon-fusion process gg —> H/A 
in which the top quark loop generates the dominant contribution. We will see that the 
production times the decay rates in the processes gg -» H/A -> tt are indeed substantial in 
a large part of the MSSM parameter space. We perform a naive estimate of the sensitivity 
that can be achieved in the search for tt resonances, a sensitivity that could allow to probe 
a significant part of the low tan (3 region of the MSSM, complementing the searches for 
t + t~ resonances that are instead sensitive to the high tan/3 region. 

The main message of the present paper is that combining the searches for heavy reso¬ 
nances decaying into t + t~ and tt pairs, and including the H —>• WW, ZZ, hh and A hZ 
channels at Ma < 350 GeV, one can possibly probe the entire [tan/3, Ma] MSSM plane 
(and not only the high-tan (3 region as is presently the case) up to large pseudoscalar Higgs 
masses; Ma values close to 1 TeV could be reached for any tan j3 with ss 3000 fb^ 1 data at 
the LHC with yfs = 14 TeV. This can be done in a model independent way and without 
relying on any additional theoretical assumption or indirect experimental constraint. 

The paper is organized as follows. We first summarise our simple parametrisation of the 
MSSM Higgs sector, further discussing and refining the /iMSSM approach. In section 3, we 
discuss the heavier Higgs production and decay rates focusing on low tan /3 and summarise 
the possible impact of superparticles. In section 4, the probing of the [tan (3, Ma] MSSM 
parameter space is discussed when all the search channels, including a projection for the 
the H/A—^ti channel, are combined first at yfs = 7+8 TeV with 25 fb^ 1 data and then at 
yfs = 14 TeV and higher luminosities. A brief conclusion is given in a final section. 
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2. The hMSSM approach 


In this section, we briefly describe the hMSSM introduced in Ref. [10], an approach that 
allows to parameterize the CP-conserving MSSM Higgs sector in a simple (as only two 
inputs are needed) and “model independent” (in the sense that we do not consider or fix 
any other SUSY parameter) way, using the information that the lightest h boson has a mass 
My l k, 125 GeV. The approach is based on several assumptions that we first summarize. 

The first basic assumption of the hMSSM is that in the basis (H ( i, H u ) of the two 
MSSM Higgs doublet fields that break electroweak symmetry, the CP-even h and H mass 
matrix can be simply written in terms of the Z and A boson masses and the angle /3 as 

M 2 _ ( M z cos2 P + m a sin2 P ~( M z + m a) sin P cos P \ ,( A M\ x A Mf 2 \ , 

$ y — (Mf + Ml) sin (3 cos / 3 Mf sin 2 (3 + M\ cos 2 (3 J y A_Mf 2 AA4 \ 2 J 


in which the radiative corrections are introduced through a 2 x 2 general matrix AA4 2 j. 
This is the usual starting point of the analyses of the neutral MSSM Higgs sector [4] and the 
calculation of the Higgs masses and couplings including radiative corrections and in which 
the SUSY scale, taken to be the geometric average of the two stop masses, Ms = , 

can be as high as a few TeV. However, if Ms is orders of magnitude higher than the TeV 
scale, the evolution from this high scale down to the electroweak scale might mix the quartic 
couplings of the MSSM Higgs sector in a non trivial way, such that the structure of the 
mass matrix at the low energy scale is different from the one given in eq. ( 2 . 1 ). 

In the ZiMSSM, we assume that the form of the CP-even Higgs mass matrix is as given 
above even at the very high SUSY scales, Ms > 0(100 TeV) that, as it will be seen later, 
are needed to consider tan/3 values close to unity 2 . 

A second basic assumption of the /iMSSM is that in the 2x2 matrix above for the 
radiative corrections, only the AA4 22 entry is relevant, AA4 2 2 > AM 2 n ,AMj 2 . In this 
case, one can simply trade AA4 22 for the by now known mass value M^ = 125 GeV using 


AM 2 M 2 (M 2 + Mf - Ml) - MjMf cos 2 2/3 
22 Mf cos 2 f3 + M\ sin 2 /3 - M 2 


( 2 . 2 ) 


This assumption is valid in most cases as the by far dominant radiative correction from 
the stop-top sector that is quartic in the top quark mass, enters only in this entry [7, 8 ]: 


AM 2 2 ~ 


3v 2 sin 2 j3 4 
8tt 2 


. Mf 
log xzf + 
mf 


X t A t 

Mf 



X t A t \ 
12 Mfy 


(2.3) 


which depends on, besides M 5 , the stop mixing parameter given by Xt = A t —fi/ tan (3 with 
At the stop trilinear coupling and /i the higgsino mass parameter. Af = \/2rnt/(v sin /3) is 
the top Yukawa coupling with v the standard vacuum expectation value v ~ 246 GeV, and 
frit the running MS top quark mass to account for the leading two-loop corrections in a 
renormalisation group (RG) improved approach. 

2 The validity of this approximation is currently studied by the LHC Higgs cross section working group 
[29]. Preliminary results in an effective two Higgs doublet model that is RG improved to resum the large 
logarithms involving Ms, have been given in Ref. [30, 31] and a more refined analysis is under way [29, 32], 
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(2.4) 


The maximal value of the h mass, M£ iax is given in the approximation above by 

Ml -> Mf cos 2 2/3 + A. M 2 22 


and is obtained for the following choice of parameters [ 20 ]: a decoupling regime with a 
heavy pseudoscalar A boson, large enough tan /3 values that allow to maximize the tree- 
level term Mf cos 2 2/3 —> M heavy stop squarks with a sufficiently large value to 
enhance the logarithmic correction log(Mf/m 2 ) and, finally, a stop mixing parameter such 
that Xf = the so-called maximal mixing scenario that maximizes the stop loops 

and hence M If the SUSY parameters are optimized as above, the maximal Mf l value 
can reach the level of M™ ax ~ 130 GeV for Ms of the order of the TeV scale, a range that 
is in general assumed in order to avoid a too large fine-tuning in the model. However, if 
tan/3 is small, the tree-level contribution Mf cos 2 2/3 to the h mass squared becomes small 
as | cos 2/3] ->0, thus requiring a substantial correction AA4 22 to obtain a sufficiently large 
Mh■ To achieve this, eq. (2.3) shows that one has to substantially increase Ms- 

In Ref. [10], the approximation AA4f 1 , A -M-12 ^ AA 1 2 2 has been checked in various 
scenarios and found to be rather good if Ms much larger than the other soft-SUSY breaking 
parameters that enter the subleading radiative corrections, such as the higgsino mass /r 
and the sbottom trilinear coupling Af, or more generally the sbottom mixing parameter 
Xf, = Ab~ /rtan/3. This assumption should be particularly justified at low and moderate 
tan/3 values where first, the bottom-Yukawa coupling A b = v2m&/(ucos/3) is not strongly 
enhanced. In the approach of Ref. [5] to parameterize the correction matrix of eq. (2.1) 
including the dominant corrections from the stop and sbottom sectors (and which has been 
used in Ref. [10] to check this second /iMSSM assumption), the entries AM.\-^ and AAij 2 
of the mass matrix when A b is set to zero are simply given at lowest order by [4, 8 ] 


A M 2 n ~ 


v 2 sin 2 /3 4 h 2 Y 2 

32vr 2 * Mg 


AMf 2 ~ 


v 2 sm 2 /3 4 
- 5, —A? x 

32vr 2 1 


jp(<i~ X,A,/Ml ) (2.5) 


They are proportional to / 1 /M 5 and hence, are small if \fi\ < Ms- Note that from the 
expressions above one can see that the two entries AAA^ and A Xif 2 are small not only for 
MsX |/r|, but also when stop mixing is small, Mg 3 > Xt- For moderate tan/3 (and also at 
large tan /3 if the sbottom corrections can still be neglected), one has At ~ Xt and the off- 
diagonal entry is further suppressed for maximal Xt = V6 Ms- Thus, the approximation 
of retaining only the entry AAA 22 for the radiative corrections should be good at least at 
low tan/3 where a very high SUSY scale is required to obtain a heavy enough h state, 
suggesting that one naturally has Ms X |/r| and eventually also Ms 2 > X t . 

In this /iMSSM approach the mass of the neutral CP even H particle and the mixing 
angle a that diagonalises the h , H states, will be given by the extremely simple expressions 

(M 2 + Mf - M 2 )(Mf cos 2 /3 + M 2 sin 2 /3) - M\M 2 cos 2 2/3 


M 2 h = 


a = — arctan 


Mf cos 2 f3 + Mf sin 2 /3 - M 2 
(Mf + M \) cos /3 sin /3 
Mfco^+Mf^Fp^M 2 


( 2 . 6 ) 


in terms of the inputs M 4 , tan/3 and the mass of the lightest h state M/j = 125 GeV. 
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(2.7) 


The mass of the charged Higgs boson is simply given by the tree-level relation 

M^± = M\ + Mr 


as the SUSY radiative corrections in this particular case are known to be very small in 
general. According to Ref. [33] where a detailed analysis of the radiative corrections has 
been recently performed, the leading one-loop correction to M^ + reads when expanding 
in powers of the SUSY scale as it is justified at low tan/3 


A MU = -YU- 9 Y_ + 0(1/Me) 

167r sin 2 0\y M^ sin 2 /3 Mj 


3a 


4 

m t 




( 2 . 8 ) 


and is therefore very small for Ms \g\- In fact, even for Ms ~ |/x| one obtains AM^± ~ 
—10 3 (250) GeV 2 for tan/3 « l(tan/3 1) and, hence, a relative correction \AM H ±/M H ±\ 
that is only about 5% (1%) for M H ± ~ 100 GeV and negligibly small for higher H± masses. 
Hence, retaining only the tree-level relation eq. (2.7) as done in the /iMSSM should be a 
very good approximation in this case. 

A third assumption of the /iMSSM is that all couplings of the Higgs particles to 
fermions and gauge bosons are given in terms of tan /3 and the mixing angle a only and, 
hence, the entire phenomenology of the Higgs particles is determined when the two inputs 
tan/3 and Ma are fixed. This means that possible corrections not incorporated in the 
mixing angle a, such as direct vertex corrections, are assumed to have a small impact 3 . In 
particular, the couplings of the neutral Higgs bosons, collectively denoted by 4>, to up and 
down-type fermions and to massive gauge bosons (including the coupling of two Higgs and 
one gauge bosons) when normalized to the SM-Higgs couplings, are simply given by: 


^ 9&u.u 

h cos a/ sin/3 

H sin a / sin /3 

A cot/3 


9<wv 

— sin a/ cos /3 sin(/3 — a) 
cos a/ cos /3 cos(/3 — a) 
tan /3 0 


g$Az/g<s>H+w- 
oc cos(/3 — a) 
oc sin(/3 — a) 
oc 0/1 


(2.9) 


The trilinear self-couplings among the Higgs bosons are also given in terms of /3 and 
a. This is clearly the case at tree-level but, to a good approximation, it remains true 
when radiative corrections are incorporated. Indeed, besides the corrections that affect 
the angle a as discussed above, the trilinear couplings receive direct corrections whose 
dominant component turns out to be simply the one that appears in the correction matrix 
AMI 2 and hence, the correction AM\ 2 °f ec l- (2-2) [34], Thus, the trilinear MSSM Higgs 
couplings are also fixed in terms of Ma, tan /3 and Mh to a good approximation. In units of 
Ao = —iM^/v, the Hhh and hhh self-couplings, which are the only ones that will matter 
for LHC phenomenology, will be then given by 

AAdooCosa 9 

A hhh = 3 cos 2a sm(/3 + a) + 3- k -- cos" a 

My sin p „ 

z 1 AM 2 


A Hhh = 2 sin 2a sin(/3 + a) — cos 2a cos(/3 + a) + 3 


oo Sill OL 9 , . 

cos 2 a (2.10) 


Ml sin /3 


3 The direct corrections can be particularly important at high tan /3 in the bottom-quark Yukawa cou¬ 
pling, as is the case of the so-called At correction which can be large if |/r| tan /3 Ms- We will show in the 
next section that even in this case, the impact of this direct correction is limited in LHC phenomenology. 
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We note that at least for the hhh self-coupling, one should incorporate the radiative cor¬ 
rections in the same approximation that has been used for the Higgs masses. This would 
be the only way to achieve a consistent decoupling limit and to make that the \hhh self¬ 
coupling indeed reaches the SM value in this limit, = 3M|/M| for a = /3 — ^. For the 
sake of consistency, one should include the radiative corrections to the other self-couplings 
in the same approximation as for A hhh- This then fully justifies the choice that we adopt 
in the /iMSSM and the expression of eq. (2.10) for the Higgs self-couplings. 

From the discussion above, one can conclude that the /iMSSM approach has two very 
interesting aspects: its economy as only two input parameters are needed to describe the 
entire MSSM Higgs sector and its simplicity, as the Higgs masses and couplings are given 
by the very simple relations eqs. (2.6)-(2.10). This would allow to considerably simplify 
phenomenological analyses of the MSSM Higgs sector which, because of the large number 
of SUSY parameters to be taken into account, rely up to now either on large scans of 
the parameter space or resort to benchmark scenarios in which most of these parameters 
are fixed. Nevertheless, the most interesting aspect of the /iMSSM is that it easily allows 
to describe scenarios with large values of the SUSY scale, Mg 1 TeV, but weak-scale 
masses for the extended Higgs sector. 

Because of the large \og(Mg/mh ) that occur, the high SUSY scale scenarios are no¬ 
toriously difficult to analyze and, before resuming the large logarithms, the MSSM Higgs 
spectrum could not be calculated in a reliable way. Until very recently, this was the case of 
the numerical tools that deal with the MSSM, such as the renormalisation group program 
Suspect [35] or the program FeynHiggs [36] that is more specialized on the MSSM Higgs 
sector, which were not reliable at too high Mg outside the decoupling regime. A new 
version of FeynHiggs in which some partial resummation of the large logarithmic terms 
is performed has become available and allows to address low tan/3 values in a somewhat 
reliable way 4 (albeit with Mh values still below Mh ~ 125 GeV at tan/3 < 2). The /iMSSM 
approach is currently being implemented in an updated version of the program Suspect. 

An immediate advantage of the /iMSSM is that it re-opens the possibility of studying 
the MSSM low tan/3 region [12], which was for a long time overlooked. Indeed, as only 
SUSY scales of the order of Mg ~ 1 TeV were assumed in the analyses performed in the 
past, one always had a too light h boson with a mass below the limit Mh >114 GeV derived 
from the negative searches of a SM-like h boson at the LEP2 collider [18,19]. For a scale 
Mg = 1 TeV, the possibilities tan /3 < 3 and tan /3 < 10 were excluded for, respectively the 
maximal-mixing scenario Xt = \/6 Mg and the no-mixing scenario Xt = 0. The situation 
became worse with the observation of the h state at the LHC and the determination that 
its mass is Mh ~ 125 GeV, i.e. well beyond the LEP limit. In fact, for MgK, 1 TeV, this 

4 This new version of the program FeynHiggs has been used to perform the following comparison [29]: for 
a (tan /3,Ma) set and given other MSSM inputs (those of the M™ ax scenario for instance), one calculates 
a , Mh and Mh using FeynHiggs and with the specific value obtained for Mh (which is in general Mh 7 ^ 125 
GeV), one recalculates the /iMSSM values of a and Mh- The relative difference between FeynHiggs and 
the /iMSSM was found to be extremely small in the entire [tan /3, Ma ] plane. Even at very low Ma and 
tan /3 values, the outputs for Mh and a differ by less than 1 %. This proves once more that the second 
assumption of the /iMSSM, i.e. that one can consider only the AM 22 radiative correction, is fully justified. 
We thank Pietro Slavich, who originally performed this comparison, for a discussion on this issue. 
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M a > M z 



M s [GeV] 

Figure 1: Contour plots in the [tan (3,M$\ plane where the there values M /, = 120,125 and 130 
GeV are obtained; the decoupling limit and maximal stop mixing are assumed. The SM inputs are 
mt = 173 GeV and as = 0.118 for the top quark mass and strong coupling constant [18]. 

relatively large My, value cannot be reached in a large part of the [tan /3, M 4 ] parameter 
space that is being explored at the LHC in the search for the additional Higgs bosons. 

Nonetheless, fixing the SUSY scale at Mg ~ 1 TeV is a very strong theoretical as¬ 
sumption and is currently challenged not only by the measured M^ [38] but also by direct 
sparticles searches at the LHC [39], especially in constrained MSSM scenarios. In the search 
for the MSSM Higgs bosons at the LHC, one would like to avoid any such assumption and 
interpret the experimental results, for instance imposing the relevant experimental con¬ 
straints in the absence of any evidence, in a way that is as model-independent as possible. 
The hMSSM approach, as no assumption on the SUSY scale nor on any other SUSY pa¬ 
rameter (except eventually that they should be smaller than Mg) is made, is more suitable 
in this respect. In fact, one is considering simply in this case a two-Higgs doublet model 
of type II [37] where the MSSM relations eqs. (2.6-2.7) are enforced; the superparticles are 
assumed to be too heavy to have an impact on the Higgs sector (as it will be shown to be 
generally the case in the next section). The only price to pay is that when the very low 
tan (3 region is addressed, one is implicitly considering a very large SUSY-breaking scale, 
making the MSSM a very unnatural and fine-tuned scenario. 

To illustrate this feature, we display in Fig. 1 contours in the [tan (3, Mg] plane in 
which one obtains the value Mh = 125 GeV for the h mass, as well as Mh = 120 and 130 
GeV. The latter examples are when one assumes that a possible mass shift of AM/, = 5 
GeV is missing from unaccounted for subleading corrections (e.g. the contributions of the 
charginos and neutralinos that we have ignored here) or unknown higher order terms (a 
theoretical uncertainty of ~ 3 GeV in the determination of My, is usually assumed [ 22 ]). 
The limit M^^Mg and maximal stop mixing Xy = \fQMg are assumed. 
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The figure has been in fact obtained from an analysis of the split-SUSY scenario where 
the large logarithms have been indeed resummed [40]. As can be seen, at high tan/3, Mg 
values in the vicinity of the TeV scale can be accommodate while in the low tan /3 region, 
extremely large values of the SUSY scale Mg are necessary to obtain = 125 GeV. 

This is particularly the case for tan/3 close to unity where a value Mg ~ 400 TeV 
is required for tan/3 = 1. The situation becomes even worse for the more natural small 
mixing situation Xt <C Mg and outside the decoupling regime when the tree-level h mass 
is reduced by Higgs mixing. In both cases, huge Mg values will be needed for tan /3 « 1 to 
reach Mj l = 125 GeV. For tan ,6 — 2. the situation is less dramatic as in the configuration of 
Fig. 1, only Mg f«20 TeV is needed to reach M/,, = 125 GeV (or the target value M^k, 120 
GeV if uncertainties are included). We thus expect that our /iMSSM approach should be 
valid down to tan/3 values as low as tan/3 ~ 2, but in our analysis we will extend the 
validity of the approach to tan /3 values close to unity. 



M a [GeV] M a [GeV] 

Figure 2: The CP-even H boson mass (left) and the mixing angle a (right) in the /iMSSM as a 
function of Ma for representative values of tan/3; the h mass is fixed to the value Mh = 125 GeV. 

Let us now illustrate the values that one obtains for the two outputs of the /iMSSM, 
the CP-even H mass Mh and the mixing angle a; the charged Higgs mass is simply given 
by eq. (2.7). These are shown in Fig. 2 as a function of Ma for several representative 
tan/3 values, from unity to tan/3 = 30. One can see that at sufficiently high tan/3 values, 
tan/3 > 10, Mh becomes very close to Ma and the angle a close to /3 — \n r, as soon as 
the pseudoscalar mass becomes larger than Ma > 200 GeV. This is a reflection of the 
well known fact that the decoupling limit, in which the A and H states are degenerate in 
mass and have the same couplings to fermions and vanishing couplings to gauge bosons, 
is attained very quickly at high tan /3. Hence, the /iMSSM approach should be a good 
approximation as it describes correctly this decoupling regime. In turn, at low tan /3, the 
mass difference Mh — Ma can be large and the angle a significantly different from /3 — 
even for Ma ~ 400 GeV meaning that the decoupling limit is reached slowly in this case. 
(For Ma ~ Mf j we are close to the regime in which the /iMSSM is not valid and one gets 
Mh~> oo and a — 7r/2; this feature will be discussed shortly). 
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tan f3 



Figure 3: Contours for the mass difference M H — Ma (left) and the coupling squared cos 2 (/3 — a) 
(right) in the [tan f3, Ma] plane. The h mass is fixed to the value Mh = 125 GeV and the red area 
at low tan /? and Ma is the one where the hMSSM is ill defined. 


This statement is made more explicit in Fig. 3 where contours for the heavier Higgs 
mass difference Mh — Ma and the square of the reduced H coupling to massive gauge 
bosons, g 2 Hw = cos 2 (/3 — a ), which is a very good measure of the distance to the decoupling 
limit. At Ma ~ 500 GeV, the difference Mh — Ma is less than 1 GeV for tan /3 > 10 while 
it is about 10 GeV for tan/3 ~ 2. However, even in this case, the mass difference represents 
about 2% of the A/H masses and, hence in view of the experimental resolution, one can 
still consider that the two states A and H are degenerate in mass. 

There is one problem with the /iMSSM approach at very low tan /3, though. If tan /3 ~ 
1 , the denominator of eq. (2.2) that expresses the dominant radiative correction AM| 2 in 
terms of the measured mass Mh becomes close to M§ + M\ — 2Af^ and, at low Ma , it 
approaches zero rendering eq. (2.2) ill defined. For tan/3 ~ 1, the pseudoscalar mass range 
Ma < 160 GeV is inaccessible, but the forbidden range shrinks to Ma < Mh for tan/3 > 3. 
As a lower bound Ma Mz has been set in the model-independent searches of the MSSM 
Higgs bosons at the LEP2 collider [19], this area in which the /iMSSM is not defined is 
therefore rather small. We will show later that this forbidden area is entirely excluded by 
the present LHC MSSM Higgs searches, in particular by and A searches as in both 
cases, the constraints can be interpreted only in terms of tan /3 and Ma and hence without 
using the /iMSSM relations of eq. (2.6). 

In fact, this “theoretically forbidden” hMSSM area is also excluded by the measurement 
of the observed Higgs boson production and decay rates at the LHC. Indeed, for these low 
Ma and tan/3 values, we are very far from the decoupling limit in which the couplings of 
the h boson are close to their SM-like values, as the LHC Higgs data in the various channels 
seem to strongly indicate [6]. However, we will refrain from using this argument to exclude 
the low Ma, tan/3 possibility, as we will prefer to perform the direct Higgs searches in a 
model independent manner, without relying on any indirect constraint. 
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3. MSSM Higgs production and decays at the LHC 

We come now to the discussion of the decays and the production at the LHC of the heavier 
A, H and H± particles in the /iMSSM. We will be mostly interested in the low tan /3 region, 
but we will first summarize the main features at high and moderate tan /3. 


3.1 Neutral Higgs decays 

At high tan/I values, say tan/3 > 10, the decay pattern of the heavier neutral H/A bosons 
is extremely simple [5,41] as a result of the strong enhancement of the couplings to down- 
type quarks and charged leptons that are proportional to tan/3, not only for the A state 
but also for the H boson. Indeed, as in the decoupling limit M A 3> Mz on has a —> /3— ); 7 r, 
the Hbb and Htt couplings normalised to the SM Higgs coupling take the limit 


9Hdd = cos a/ cos /3 


Ma>M z 


tan /3 = g A dd 


(3.1) 


The neutral = A/H states will decay almost exclusively into t + t~ and bb pairs, 
with branching ratios of BR(4> —>• rr) ~ 10% and BR(<F bb) ~ 90%. This is a simple 
consequence of the fact that the partial widths are proportional to respectively ( m T tan/3 ) 2 
and 3(m;,tan/3 ) 2 when the color factor is included; rn T = 1.78 GeV and the MS bottom 
quark mass defined at the scale of the Higgs mass is nib ~ 3 GeV, implying thus, 

BR(4> —»• rr) « 1 — BR($ —> bb) m 2 /(m 2 + 3m 2 ) « 0.1. (3.2) 


At high tan/3, all other decay channels of the H/A states are strongly suppressed. 
This is particularly the case of the decays into top quark pairs, despite of the large value 
m.t^>mb, as the Higgs coupling to up-type quarks are inversely proportional to tan/3, 

9Huu = sin a/ sin /3 Ma zZ1/ Iz \j t an /3 = g Auu (3.3) 


rendering very small the = H, A partial widths, given in the Born approximation by 

T(4> -> tt) = 3G F m t 2 /(4^2vr) x g% tt /3f (3.4) 

where /3* = (1 — 4/ra 2 /Af|) 1 / 2 and p = 3 (1) for the CP-even (CP-odd) Higgs boson. 

This is also the case of Higgs decays involving gauge and Higgs particles in the final 
state. In particular, one should have in principle also the decay modes H —> VV with 
V = W, Z and H —» hh in the case of the CP-even and A —> Zh in the case of the CP-odd 
Higgs bosons. However, the partial decay widths of the H particle into massive gauge 
bosons F(H -> VV) are proportional to the square of the reduced coupling 

9hvv = cos(/3 - a) —> 0 = g A vv (3.5) 

which becomes zero in the decoupling limit as is the case for the pseudoscalar A boson, 
that has no tree-level couplings to VV states as a result of CP-invariance. For the latter 
state, the possibility A-^hZ for M A > M^+Mz~220 GeV, i.e. near the decoupling limit, 
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will have a suppressed rate as the coupling gAhZ = 9hvv tends to zero at large tan/3. 
Indeed, an expansion in terms of 1/M\ gives [12] 

9HVV = 9Ah.z Ma ^ Z Mf/(2M|) x sin4/3 - AM| 2 /(2M|) x sin 2/3 (3.6) 


and, at high tan/3, both sin4/3 and sin2/3 are proportional to cot/3 so that the limit 
9hw —> 0 is reached faster in this case. The same is true for the decay H -A- hh when 
Mm > ‘IMh as the trilinear Higgs coupling of eq. (2.10) for Mh > 2 Mh reaches the limit 

9Hhh Ma ^ Z —3AMf 2 /(2Mf) x sin 2/3 (3.7) 

and is thus strongly suppressed at high tan/3 that implies sin 2/3 = 2 tan/3/(l+tan 2 /3) —>-0. 

The situation is drastically different at low values of tan (5 when the heavy Higgs states 
are kinematically allowed to decay into top quark pairs, Mh ~ M^> 2mt [5]. Indeed, 
H/A^-tt become by far dominant g$ t t oc mt/ tan /3 is so strong that it leaves no chance to 
the other possible channels. This is clearly the case for the H/A -A bb, t + t~ rates which 
become negligibly small as the couplings g$dd are not enhanced anymore and mt 3> mb, m T . 
This is also the case of the decays A -A hZ and H -A VV at large since the couplings 
approach zero in this case. For the H -A- hh decay, there is still a component of the gHhh 
coupling of eq. (3.7), the one oc sin 2/3 ~ 1 for tan /3 ~ 1, that is non-zero in the decoupling 
limit. However, besides the fact that the Htt coupling is larger than the Hhh coupling, 
the partial decay width for the process H -A- hh decreases as 1/Mm 

T(H -a hh) = (G F Ml) /(16V27T Mh) x X 2 Hhh (l - 4 Ml/M 2 H ) l/2 (3.8) 


while T(H -A-tt) oc Mh and, hence, largely dominates at high Mh (slightly above M^~2mt 
however, BR(i7 — hh) stays significant as T(H^-tt) is suppressed by /3f near threshold). 

The situation is opposite when the decays H -A VV are considered. Indeed, because of 
the contributions of the longitudinal components of the V = W,Z bosons that grow with 
the energy scale, the partial decay widths increases as Mfj and not as Mh 


T(H -> VV) = 


GfSv 2 n,f3 a ( 1 ^My 

if* 1 """* V 14 


+ 


12 My 

~Mf 


M a»M z Gf^v 


9 2 hvvM 3 h (3.9) 


where again /3y = (1 — AMy/Mjj) 1 / 2 and 5y = 2(1) for V = W(Z). At low tan/3 and 
high Mh values, one of the components of gnw given in eq. (3.6) (the one oc sin 4/3) 
vanishes while the other component tends to gHW —^>AM 2 2 /Ma. Because of the 
enhancement of the decay rate by M\, one would have then a partial width T(H -A VV) 
that is suppressed by a power 1 /Mh only and hence, does not become completely negligible 
compared to H —> tt even at very high Mh- For instance, at Mh ~ 500 GeV, the branching 
ratios for the decays H -A- WW and H -A ZZ are still at the 2% and 1% level respectively. 
This is appreciable and, at least in the case of the ZZ decay, it is of the same order of 
magnitude as the branching fraction of the observed 125 GeV h boson, with the advantage 
that the ZZ pair has a much larger invariant mass with a significantly smaller background 
(which compensates for the smaller Higgs production cross section as will be seen later). 
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If the H/A states have masses below the 2 mt kinematical threshold, the two-body 
H/A —> it decays are not allowed anymore. Off-shell three-body decays A/H —> it* —> ibW 
are possible, but the rates are suppressed by an additional electroweak factor and the 
virtuality of one of the top quarks [42], The gauge and Higgs decays of the H/A states 
would then become significant at low tan/3 values. In the mass range M^+Mz < Af 4 < 2 mt, 
the decay A —> hZ will be dominant: the reduced coupling gAhZ = Shw A/| /M\ is 
only moderately suppressed and the full AhZ coupling is still substantial compared to the 
tiny Abb coupling. Likewise, for 2 Mh < A Ih < 2 mt, the decay mode H —> hh is dominant 
as the coupling guhh at low tan /3 will stay appreciable. The two-body decays into massive 
gauge bosons H —> WW and ZZ are also significant below the 2 mt threshold. 

The bosonic decays will also be non-negligible at intermediate values of tan /3, tan /3 ~ 
\/fnt/mb ~ 5-10, when the A/H couplings to top quarks are suppressed while those to 
bottom quarks are not yet strongly enhanced. However, below the 2 mt threshold, when 
the Higgs couplings to the bosonic states are not too suppressed and the only competition 
will be due to the T —> bb decays that is only slightly enhanced, the rates for the H —> 
WW, ZZ, hh and A —>• Zh channels will be smaller than at low tan /3 values. 

The branching fractions for the various Higgs decays discussed above are displayed in 
Fig. 4 in the [tan /3, Ma] plane assuming the LMSSM with Mh = 125 GeV. The Fortran 
program HDECAY [41] where the hMSSM relations were implemented has been used. The 
color code is such that the red area is when the considered decay rates are large, while the 
blue area is when they are small (for alternative and more easily readable figures for the 
branching ratios, as well as for the production cross sections, see for instance Refs. [10,12]). 
The white areas are when the decay rates are very small, below the minimal value of the 
scale in the color axis. As can be seen, the H/A —> tt decays are important at high tan/3 
values. The branching ratios for the decays H/A —>• bb follow that of tt final states when 
a factor 9 is included and are, hence, largely dominant. The decays H/A —>• ti are by far 
leading at low tan (3 for Ma,h ^ 350 GeV (one notices that at least for A, they are also 
significant slightly below the ti threshold). The bosonic decays H —> WW, ZZ, hh and also 
A —> hZ have reasonable rates only for tan (3 < 3 and below the 2 mt threshold. 

A final word should be devoted to the total Higgs decay widths, which are displayed 
for the A and H particles in Fig. 5, again in the [tan/3, Ma] LMSSM plane. In the low and 
high tan (3 regimes, one can consider only the dominant fermionic decays of the 4> states 
and, to a good approximation, the total decay widths read (up to phase space factors) 

T tot _ r tot(_ r tot ± ) _ G F A/$/(4^2vr) [(ml + 3m 2 b ) tan 2 /3 + 3m? cot 2 /3] (3.10) 

For tan/3 ~ 1 and tan/3 ~ 60, one obtains a total decay width that is approximately 
r^ ot ss 5%Af$ and, compared to the Higgs mass, it is not very large. Hence, to a good 
approximation the A/H states can be considered as narrow resonances in most cases. 

Note that for the H state, as the branching ratios and total decay width are shown as 
a function of Ma, some peculiar features can be observed. These are explained by the fact 
that there is a large splitting between Mh and Ma at low tan /3 and Ma values which lead 
to, for instance, the opening of the H —» ti mode already at Ma < 200 GeV and, hence, 
suppressed H —> WW, ZZ decays but a large total decay width for tan/3 ~ 1. 
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Figure 4: The branching ratios of the neutral Higgs bosons in the [tan /?, Ma] parameter space of 
the /iMSSM with the constraint Mh = 125 GeV. Shown are the rates for the decays A/H —> rr 
(top), A —> hZ and H —> hh , WW, ZZ (middle) and A/H —>• tt (bottom). 
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Figure 5: The total decay widths of the neutral A (left) and H (right) bosons in the [tan/3, Ma\ 
parameter space of the 6MSSM with the constraint Mh = 125 GeV. 


3.2 Neutral Higgs production 


Let us turn now to the production of the neutral MSSM <f> = H/A bosons at the LHC. Also 
in this case, the cross sections crucially depend on the considered tan (3 regime and in most 
cases, the two processes that play a leading role are the gluon-fusion mechanism gg —> $ 
which is initiated by a heavy quark loop [43] and the associated Higgs production with 
6 -quarks, gg/qq — > 66 +<I>, which at high energies and if no-additional 6 -quark is considered 
in the final state, is equivalent to the fusion process bb —> [44], All other processes, in 
particular vector boson fusion and associated production with a massive gauge boson for 
the CP-even H state, qq —>• Hqq and qq —>• HV, and associated production with top-quark 
pairs for both the H and A states, pp —> 4>tf, have much smaller rates as the couplings 
9hw and g$ t t are suppressed and/or the available phase space is not favorable. 

At leading order in perturbation theory, the partonic cross sections for the bb —> and 

gg —> processes can be written in terms of the partonic c.m. energy s and as [3] 

a(bb ->■ $) = ^ g\ bb 6{s - Ml) (3.11) 


Hd9 -> $) = 0 ^ 0 F( ^ Af|<5( 


(3.12) 


In the case of the gg —> process, the quarks Q running in the loop should be taken to 
be the heavy bottom and top quarks with Higgs couplings given in eq. (2.9) and masses 
incorporated into the reduced variables tq = Ml/ArriQ. The form factors for spin-^ 
fermion loops in the case of a CP-even H and a CP-odd A bosons are given by 

A i/ 2( t ) = 2 \ t + ( r - 1 )/(^ r )] r_2 

4%(r) = 2 r_ 1 /(r) (3.13) 

where the function f(r) above and below the r = 1 kinematical threshold is defined as 


f( T ) = < 


arcsm 


1 + y/l - t 1 
log - , — ITT 

1-Vl-T- 1 


1 2 


T < 1 
r > 1 


(3.14) 
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While the amplitudes are real for M< j> < 2rriQ, they develop an imaginary part above the 
kinematical threshold. At very low Higgs masses, compared to the internal quark mass, 
the amplitudes for a scalar and a pseudoscalar states reach constant but different values 

M<i> <C 2viq : A^ 2 (tq) — > 4/3 , Ai/ 2 { t q) 2 (3.15) 

Instead, in the opposite limit, M, j> 3> 2 uiq, chiral symmetry holds and the amplitudes for 
the CP-even and a CP-odd Higgs bosons are identical (as in the 66 —>■ case), 

Ms > 2 m Q : Af /2 (r Q ) -> -[log(4r Q ) - mf/{2T Q ) (3.16) 

The maximal values of the amplitudes occur slightly above the kinematical threshold where 
one has for the real parts R e(A^ 2 ) ~ 2 and Re(Ajy 2 ) ~ 5. 

At high tan/3, the strong enhancement of the Higgs couplings to 6-quarks and the 
suppression of the couplings to top quarks and gauge massive bosons makes that only these 
two processes are relevant, with the gluon-gluon fusion mechanisms dominantly generated 
by the bottom quark loop. The cross sections a(gg —> <h) and a(bb —> <f>) are of the same order 
of magnitude and can be so large that they make the process pp—> gg+bb —>■<!>— >t + t~ , with 
the branching fraction for the decay 4> — >t + t~ being of order 10% as seen previously, the 
most powerful LHC search channel for the heavier MSSM Higgs bosons. The pp —> bb+H/A 
mode with A/H —> 66, which has an order of magnitude larger rates in principle (if no 
high-pr 6-quark from production is required), has also been considered but the sensitivity 
is smaller as this fully hadronic process is subject to a much larger QCD background. 

At high tan /3, the production rates are approximately the same for the H and A states 
in both the 66 and gg fusion cases as discussed earlier. While a(bb —> <h) is known up to 
NNLO in QCD perturbative theory [45], cr(gg —> 4>) is instead known only up to NLO in 
the limit M<j> > 2niQ that we will be mainly interested in [46]. For the top-quark loop, we 
will nevertheless also include the NNLO QCD corrections [47] that are in principle only 
valid for < 2 viq as advocated in Ref. [48]. The precise values of the cross section times 
branching fractions cr(pp —>• 4?) x BR(4>— >t + t~) for a given [tan/3, Ma\ MSSM point have 
been updated in Refs. [48,49] and the associated theoretical uncertainties from missing 
higher order perturbative contributions, the parametrisation of the parton distribution 
functions and uncertainties on the inputs a s and 6-quark mass, have been estimated to be 
of the order of 25%. Any effect below this level, such as the SUSY effects that we will be 
discussed later in this section, should be considered as small. 

Again, at low tan/3, the situation is very different. The cross sections for the 66 —> $ 
processes are now very small as the <F66 coupling is not enhanced anymore. For Ma <1 200 
GeV, this is also the case of the associated production with tt pairs as a result of a the 
small phase-space and, in the case of the H state, of the vector boson fusion qq Hqq 
and associated production with a gauge boson qq —> HV as a result of the suppressed 
HVV coupling (the A state cannot be produced in these two processes as there is no 
AVV coupling). The only process which would have a reasonable production cross section 
would be the gluon-fusion process gg —> with, this time, the leading contribution being 
generated by loops of top quarks that have couplings to the Higgs bosons that are only 
slightly suppressed compared to the SM Higgs case. 
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Figure 6: The production cross sections of the MSSM heavier Higgs bosons A (left) and H (right) 
at the LHC with y/s = 8 TeV (top) and 14 TeV (bottom) in the [tan 0. Ma] ZiMSSM plane; the 
channels gg —>■ H/A and bb —> H/A have been considered. 
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The production cross sections <r(gg —>• 4 - X ) and a(bb —> $ + X ) with $ = A (left) 
and Q = H (right) are displayed at the LHC in the [tan/3, Ma] /iMSSM parameter space 
for y/s = 8 TeV (top) and 14 TeV (bottom); the MSTW parton distribution functions [50] 
have been used. The rates for bb —> <h have been obtained using the program SusHi [51] 
and an adapted version of HIGLU [52] has been used for gg —> <h. As can be seen, the cross 
sections are rather large in particular at high tan/3 and, for gg —> 4>, also at low tan/3 
when the Higgs couplings to b~ or t-quarks are strong and at relatively low Ma when the 
phase space is not too penalizing. Even for Ma = 500 GeV (1 TeV), the production rates 
are significant at y/s = 8 TeV (14 TeV), if tan/3 is sufficiently high or low. 

3.3 The case of the charged Higgs boson 

A final word should be devoted to the case of the charged Higgs boson, whose coupling to 
fermions is proportional to 

9H±ud oc m d tan/3(l + 75 ) + m n cot/3(l - 75 ) (3.17) 

The coupling is large at low tan/3 when the component rn u j tan/3 is not suppressed and at 
very large tan (5 when the component tan /3 is enhanced, so that many aspects discussed 
for the pseudoscalar Higgs boson hold also in this case [5]. At low mass, M H ± < 160 
GeV, which corresponds to Ma < 140 GeV, the charged Higgs boson can be produced in 
the decay of top quarks that are copiously produced at the LHC, gg + qq —> tt with one 
top quark decaying into the dominant t —)• bW mode and the other into t — > bH + . For 
M h ± ~ 140 GeV, the latter channel has a branching ratio ranging from order ~ 10% for 
tan/3 ~ 1 or tan/3 ~ 60 to order ~ 1% for tan/3 ~ 7-8 when the Higgs couplings are the 
smallest. In this low mass range above, the H± boson will decay almost exclusively into 
tv final states but some some competition with the hadronic decay channel H + —> cs will 
occur at very low tan /3. The branching fractions are shown in Fig. 7 as a function of 
M h ± for two representative tan/3 values, tan/3 = 2 and tan/3 = 30. 
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Figure 7: The branching ratios of the charged Higgs boson as a function of its mass for tan/3 = 2 
(left) and tan/3 = 30 (right); the hMSSM with the constraint Mh = 125 GeV is assumed. 
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o(pp -> H*t(b)) (pb) 


At higher masses, the state will be mainly produced in the three-body production 
process pp — > tbH ± which, at high energies, is equivalent to the two-body channel gb —> 
H±t if no additional final state b quarks are detected [53]. Again, significant rates occur 
only at very low or very large values of tan/3 when the H±tb coupling of eq. (3.17) is large 
(some small additional contributions from the tree-level qq —>• H + H~ and loop induced 
gg —y H + H~ pair and associated qq^H^+A/h/H production modes are also possible 
[5]). The cross sections have been derived in Ref. [54] where the two possibilities for the 
process, pp —> tbH~ and gb —)• tH ~, are properly matched and some numerical grids have 
been provided for the MSSM. The output of these grids for the production rates at = 8 
and 14 TeV is shown in Fig. 8 in the [tan (3,M H ±\ plane. 

At high tan /3, as shown in Fig. 7, the decay branching fractions are BR (H + — >• 
tv) ~ 10% and BR (H + — > tb ) ~ 90% exactly for the same reasons discussed previously 
for the H/A particles. All the other decay channels can be safely ignored so that the 
main search channel would be pp —> H±t{b) production with H + —» tv. As will be seen 
later, the process is however less powerful in probing the MSSM parameter space than the 
pp —?• H/A —>• rr channel discussed earlier. 



Hf(GeV) 



Mt(GeV) 


Figure 8: The matched production cross sections pp — > H + t(b) at the LHC with ^fs = 8 TeV 
(left) and 14 TeV (right) as a function of M H ± for several tan/3 values; the numbers were taken 
from the grids given in Ref. [54]. 

At low tan/3, the dominant decay channel will be by far the H + —>• tb mode with a 
branching ratio close to unity for a sufficiently heavy state when phase-space effects are 
irrelevant (one should note though, that slightly below the mt + m & kinematical threshold, 
the three-body decay with an off-shell top quark, H + —> bt* —> bbW , is also important in 
the case tan/3 ~ 2 as shown in Fig. 7). The main search channel in this case would be the 
pp — > i7 ± t(6) —> ttbb mode. This process, which is sensitive to the same area of [tan /3, Ma] 
parameter space as the H/A —> tt channels discussed before i.e. the low tan/3 and high 
Ma regions, has been considered in the past and found to be of limited use as it is subject 
to a large QCD background [55]. However, a recent CMS analysis [56] gave interesting and 
more optimistic results that we will discuss in the next section. 
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Finally, at masses above M H ± > 160 GeV, an interesting decay channel would 
occur, namely H± —tWh. Nevertheless, and in contrast to the similar A —> hZ decay mode 
discussed previously, this channel has to compete from the beginning with the dominant 
—y tb decay. Only at moderate tan /3 and low M H ± that this decay has a sizeable 
branching ratios, of the order of a few percent, allowing for H± searches in the interesting 
channel pp —> tbH ^ —> tbWh which, experimentally, has not been considered so far. 


3.4 Impact of the SUSY spectrum and dark matter 

An important question would be if the MSSM Higgs production times decay rates are not 
affected by the presence of supersymmetric particles. These could have two impacts: first, 
they could contribute virtually to the processes and modify the production cross sections 
and decay branching ratios. This issue is directly related to the third assumption of the 
/iMSSM, namely that the couplings of the Higgs bosons are simply given by eqs. (2.9)- 
(2.10) and no direct correction is involved. Second, SUSY particles could appear in the 
decays of the Higgs particles and alter the branching ratios for the standard channels that 
are searched for. This possibility would also invalidate the simple /iMSSM approach as 
some SUSY parameters would be then required to describe Higgs phenomenology. Both 
issues have been discussed e.g. in Ref. [12] and below, we simply summarise the main 
points with details concerning the low tan/3 region and the decays H/A —> tt. 

For what concerns the production processes, besides the standard top and bottom- 
quark loops, there are also squark (and mainly stop) loops [57] that contribute to the 
production of the CP-even H boson in the gluon-fusion channel, gg^-H] the CP-odd A 
states does not couple to identical sfermions and there is no-squark contribution to gg —> A 
at lowest order. However, as the Higgs-squark couplings are not proportional to squark 
masses, the contributions are damped by powers of the squark mass squared oc 1 /tHq and 
should be small for sufficiently heavy squarks. This is particularly the case at high tan /3 
values where the standard bottom-quark contributions are so strongly enhanced that the 
impact of squarks becomes negligible. At low tan/3 values, as one needs a large SUSY scale 
Mg 2 > 1 TeV in order to accommodate an h boson with a mass Mh ~ 125 GeV, the impact 
of the too heavy squarks should also be negligible in the gluon-fusion process. Hence, in 
most cases, these SUSY loop contributions can be ignored in the production modes. 

SUSY particles can have a large impact also through the Higgs boson couplings. In¬ 
deed, besides the radiative corrections that affect the Higgs mass matrix eq. (2.1), there 
are additional one-loop vertex corrections that modify the Higgs-fermion couplings and 
which are not described by the /iMSSM. These corrections are in general only important in 
the case of 6-quarks and only at high-tan (5 and large p values, since they grow as /rtan/3. 
The dominant components are due to the contributions to the Higgs-66 vertices from the 
strongly interacting sbottoms and gluinos and the weakly interacting higgsinos with top 
squarks. They can be approximated by [58] 


A h ^ 


2 a 


A? 


—mg/max(mg,m- b _) + A t /max(p , m~ 


/rtan/3 


(3.18) 
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They affect mainly the heavier Higgs couplings that become in the limit Ma Mz 


3Hbb ~ 9Abb ~ 9H±tb 


tan/3 

1 + Aft 


at high tan/3 


(3.19) 


For the lighter h state, the coupling ghbb is not affected in this limit and stays SM-like. 

Nevertheless, as already discussed in many places including Refs. [12,49], this correc¬ 
tion has only a limited impact in the case of the full pp -A = H/A -A rr process as the 
correction appears in both the production cross section 


cr(gg + bb -> <f>) oc (1 + A b ) 2 


and in the rr decay branching fraction, 

BR(4> -> rr) = T(4* -> rr)/[(l + Aft)" 2 r($ -> bb) + T(<h -> rr)] 

, and it largely cancels out in the product of the two 

a x BR ~ 1 — -Aft (3.20) 

5 

Hence, only when the Aft correction is huge and larger than unity (a feature that might 
put in danger the perturbative series) that its impact on the pp -A rr cross section times 
decay rate is of the order of the theoretical uncertainty, about 25% as discussed earlier. At 
low tan/3 values, and eventually also at intermediate tan/3 values, the Aft correction is not 
enhanced and its effects should be rather small. Nevertheless, this is not the case of all 
process and, in particular, the search channel pp—^H/A with H/A^bb in which the Aft 
impact is in fact doubled in the production times decay rates, and the pp -A H^bt mode 
with H + -A tv at low M H ± and H + — > tb at high M H ±. We will see, however, that these 
processes do not play a leading role in MSSM Higgs searches at the LHC at high tan /3. In 
summary, and to first approximation, one can thus consider that the Aft correction has a 
limited impact on the Higgs searches in the /j.MSSM. 

For the second option, namely that light SUSY particles could contribute to the decays 
of the Higgs bosons, the situation is also relatively simple; see Ref. [5] for a review. At 
very high tan /3, the partial widths of the H/A-+ bb, t + t~ as well as H + -» tb, t + v decay 
modes are so strongly enhanced, that they leave no room for the SUSY decay channels. 

At low tan /3, high values of the SUSY scale are required, resulting in large squark 
masses (at least in universal models in which the squark masses of the three generations 
are related) that make the Higgs decays into squarks kinematically closed for reasonable 
Mj 4 values. If the masses of the sleptons are disconnected from the SUSY scale and are 
made small enough for the decays of the heavier Higgs bosons into slepton pairs, H —> £i£j 
and A -A (again, CP invariance forces the A boson not to couple to identical sfermions), 
to occur. Nevertheless, except for the Higgs-stau couplings at sufficiently high tan /3 values 
when the competition from the standard channels is though, the Higgs-sleptons couplings 
are in general small making these channels very rare and their impact limited. 
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Thus, only decays into charginos and neutralinos could play a role and affect signifi¬ 
cantly the Higgs branching fractions in the standard channels' 5 . Let us briefly comment on 
these channels. 

Three conditions must be fulfilled in order to have significant rates for Higgs decays 
into charginos and neutralinos, H/A —> XiXj (with i,j = 1 - - -4), H/A —> xfxj (with 
i, j = 1,2) and H ± -A X fxj [5,60]. 

i ) One needs that some of the x states are light, > 2 m x , in order to allow for some 
decay channels to be kinematically open. 

ii) One needs to have significant &XX couplings; these couplings are maximal when 
the x final states are mixtures of higgsinos and gauginos, a feature which requires 
comparable higgsino and gaugino mass parameters, [i « M 2 . 

Hi) One needs that the standard Higgs decay modes are not enhanced and hence, not 
too low or too large values of tan /3 where, respectively, the Higgs-top and the Higgs- 
bottorn couplings are enhanced. 

The maximal Higgs decay rates into charginos and neutralinos are obtained at moder¬ 
ate tan f3 when all X X channels are kinematically accessible. In this case, as a consequence 
of the unitarity of the diagonalizing chargino and neutralino mixing matrices, the sum of 
the partial widths do not involve any of the elements of these matrices in the asymptotic 
regime M, j> 2 m x where phase space effects can be neglected. The sum of the branching 
fractions of the three Higgs bosons = H/A/H± decaying into the various XX final states 
is then simply given by {6w is the electroweak mixing angle) [60] 

\ = _ (l + | tarn 6 W ) Mjy _ 

' * 1 (l + 5 tan 2 @w) Myy + rhf cot 2 /3 + (m 2 + m 2 ) tan 2 f3 

when only the leading tt, bb and rr modes for the neutral and the tb and tv modes for 
the charged Higgs bosons are included in the total widths. This is approximately the case 
as we are close to the decoupling limit when these SUSY channels are accessible and the 
other standard decay modes such as H — > VV, hh and A — > hZ are suppressed. 

The branching ratios when all ino states are summed up are shown for the three MSSM 
Higgs states in Fig. 9 as a function of tan/3 for Ma = 600 GeV. The hMSSM relations for 
the Higgs sector have been enforced and the other relevant SUSY parameters are fixed to 
// = M 2 = 200 GeV, assuming that the bino and wino soft SUSY-breaking parameters are 

5 Here we will consider models in which the neutralino yj is the lightest SUSY particle (LSP). In gauge 
mediated SUSY breaking models (GMSB), a very light gravitino can be the LSP and decays of the heavier 
Higgs bosons into a gravitino and a neutralino or a chargino, H/A —> G\i and H ± -A- Gxt> are possible 
(Higgs bosons do not couple to pairs of gravitinos). However, as discussed in Ref. [59], the partial widths 
for these decays are inversely proportional to the square of the SUSY breaking scale, Mj = (m(jMpi anc k) 2 , 
and need Ms < 0(few 100 GeV) to be substantial. This possibility is ruled out not only from direct SUSY 
particle searches but also from the fact that the large value of the MSSM h boson mass, Mh = 125 GeV, 
requires a SUSY scale in the multi TeV range in GMSB scenarios (Ms > 3 TeV in minima GMSB) [38]. 
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related by the unification condition M 2 ~ 2M\. One can see that the branching ratios for 
the three Higgs particles are indeed similar and that they do not dominate at low nor at 
high tan/3. For instance, they are less than 25% (which is the magnitude of the theoretical 
uncertainty on the production cross sections) for tan /3 < 2 and tan /3 > 30 as can be seen 
from the figure. In contrast, the XX branching rations can be large for intermediate values 
of tan/3 when the Higgs couplings to top (bottom) quarks are suppressed (not strongly 
enhanced) and, for instance, they reach the level of ~ 70% at tan /3 ~ 5-10. 



Figure 9: The sum of the branching ratios for the decays of the three heavier MSSM Higgs bosons 
into charginos and neutralinos as a function of tan /3 for Ma = 600 GeV. 

Nevertheless, this possibility with large Higgs decay rates into SUSY particles seems 
unlikely. First, low and comparable values of the wino and higgsino mass parameters 
M 2 ~ /j, < 300 GeV that would lead to light charginos and neutralinos in the decay products 
of not too heavy H/A/H ^ bosons, are constrained by the direct searches for these particles 
at LEP and the LHC [18,61,62], In particular, the associated production of the lighter 
chargino and the next-to-lightest neutralino qq —> xtx 2 wou ld lead to large cross sections 
at yG = 8 TeV and the decays xf —> Wxi and X 2 ^Xi> with leptonic gauge boson 
decays, would have significant branching fractions. The search for leptons plus missing 
energy at the first run of the LHC, in particular the clean trilepton events from the chain 
pp X1X2 WZ X °iXi imposes severe restrictions on the parameter space. 

For instance, for M\ ~ 100 GeV, the area fi « M 2 < 200 GeV that corresponds to the 
choice adopted for Fig. 9 is by now excluded by the LHC data 6 [62], 

6 Note that we have adopted the same choice of gaugino-higgsino parameters as in the benchmark scenario 
of Ref. [21] that has been used to interpret the experimental limits in the pp — > tt searches made by the 
ATLAS and CMS collaborations [16,17]. Hence, this choice leads to a large branching fraction for Higgs 
decays into \X states for Ma > 300 GeV and weakens the experimental constraints that can be obtained 
from H/A searches in the pp tt process, while it is apparently excluded by direct SUSY searches. 
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To evade these experimental bounds, one needs either to increase the parameters /U2~ 
M2 well above 200 GeV with a consequence that the phase space for the Higgs decays 
will be limited, or to make that the light x states are either pure higgsinos (/jCMj) or 
pure gauginos (|/> M 2 ), which then suppresses the Z\X couplings e.g. and hence the 
trilepton signals, leading to Higgs couplings to the kinematically accessible charginos and 
neutralinos that are too small. Thus, in all these cases, the SUSY decays are suppressed 
and do not jeopardize the Higgs signals in the standard search channels. 

These arguments become stronger in the case where the lightest neutralino Xi is forced 
to be the candidate for the dark matter in the universe, with a cosmological relic density as 
measured by the WMAP/Planck teams, 0.09 < Qh 2 < 0.12 [63] with h being the reduced 
Hubble constant. Traditionally, four regions of the MSSM parameter space have been 
advocated to fulfill this condition for the LSP neutralino [64] and we list them below. 

i ) The so-called “well tempered neutralino” region, with a mixed gaugino-higgsino 
LSP that allows for significant couplings to gauge and Higgs bosons; this leads to a good 
LSP annihilation rate into bosonic finals states, X1X1 WW, ZZ,hZ. For low M 2 and 
/x values, this is the region discussed above that is constrained by the multi-lepton plus 
missing energy searches. 

ii) The bino like neutralino region, /x Mi, where the main LSP annihilation channel 
is into third generation tau leptons, with the exchange of light sleptons in the t-channel, 
X1X1 t + t~ , and the f co-annihilation region in which the lightest f\ state is almost 
mass degenerate with the LSP, m^ ~ m x 0, and the correct relic density is provided by the 
process fifi —> SM particles. If the SUSY scale is high and a kind of universality is assumed 
for sfermions, the fi state will be too heavy and both channels become inoperative. 

in) The regions where the LSP is almost a pure higgsino or gaugino state and hence 
with small couplings to the H/A/H± bosons. The correct cosmological relic density is 
provided by the co-annihilation of the mass degenerate xf and X2 states that need to be 
very heavy (above 1 TeV) and inaccessible in the decays of TeV Higgs bosons. 

iv) Finally, there is the Higgs-pole region [65] in which an efficient LSP annihilation 
into SM particles is provided by the exchange of an almost on-shell A boson in the s— 
channel; one thus needs M4 ~ 2m x o (the possibility of /x-boson exchange [66] leading to 
m x o ~ 60 GeV is by now unlikely). In the past, the high tan (3 region was favored and the 
most discussed annihilation channel was X 1 X 1 —> A —> bb. At low tan (3, a new possibility 
opens up, the channel X 1 X 1 —> A —^ tt which can also lead to the correct relic density. 

In Fig. 10, we display the areas of the [M2,/x] parameter space in which the relic 
density of the lightest neutralino, calculated using the program micrOMEGAs [67], is as 
determined by the WMAP/Planck collaborations, i.e. 0.09 < Qh 2 < 0.12 [63]. We have 
chosen M4 = 500 GeV and tan (3 = 2 for the LMSSM inputs, and assumed very heavy 
sfermions. In the left plot we have fixed the bino mass to M\ = 100 GeV, while in the 
right plot, we used the unification condition M2 = 2 M\. Outside the areas excluded by 
LEP and LHC ino searches at Mi = 100 GeV ~m v 0, only two areas lead to a correct relic 
density: small M2 or /x values ( < 150 GeV) that allow for X1X2 an d X1X1 co-annihilation 
or a mixed bino-higgsino LSP. Another area opens up if m x 0 ~ |M2 ~ \Ma a s can be 
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seen in the right-figure: the A funnel in which the LSP efficiently annihilates through 
the channel X1X1 —> A —>• tt. In most of this area, the decay A —> X1X1 is kinematically 
closed (or phase-space suppressed) and, because Xi is the LSP, so are all Higgs decays into 
superparticles. 

Hence, in all cases, the requirement that the lightest neutralino is the dark matter in 
the universe with the correct relic density makes that the decays of the Higgs bosons into 
charginos and neutralinos should not occur, or at least should not dominate. 



Figure 10: Points in the [M 2 , fi\ plane which satisfy the dark matter constraint with Ma = 500 GeV 
and tan /3 = 2 for Mi = 100 GeV (left) and with the unification condition for the gaugino masses 
Mi = O. 5 M 2 (right). In the left figure, the parameter space excluded by the LEP and LHC direct 
searches of charginos and neutralinos is indicated. 

One concludes from all the discussions of this subsection that it is rather unlikely that 
the SUSY particles make a significant impact in the phenomenology of the MSSM Higgs 
bosons, either in their virtual contributions to the production and/or decay processes (in 
particular since the SUSY effects should be larger than the ~ 25% theoretical uncertainty 
that affects the production rates) or in the direct appearance in the decays of the heavier 
Higgs states. One can thus assume that the superparticles are very heavy and/or too 
weakly coupled and that they decouple from the MSSM Higgs sector, except of course in 
the radiative corrections to the CP-even mass matrix eqs. (2.1). 

This is, in fact, another way of stating the third assumption of the /iMSSM discussed 
in section 2 and one can thus consider this effective approach as a very good benchmark. 
Ignoring the SUSY effects is a rather reasonable attitude since, besides the tremendous 
simplifications that the /iMSSM introduces in the description of the Higgs sector, it leads 
to a straightforward interpretation of the experimental constraints, that do not need to be 
”de-convoluted” from these complicated effects when they are included. 

Nevertheless, the possibility of light charginos and neutralinos with sizable couplings to 
the Higgs bosons cannot be totally excluded at the moment. On should therefore continue 
performing searches for the H/A/H + XX SUSY decays and other more direct searches for 
these particles such as the tri-lepton signal events. In the full MSSM with light SUSY 
particles, all these processes provide complementary information. 
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4. The probing of the MSSM parameter space 

4.1 Interpretation of the fermionic Higgs decay modes in the hMSSM 

As discussed earlier, the most efficient channels that allow to probe the MSSM parameter 
space at the LHC are the search for charged Higgs bosons coming from top quark decays, the 
process t —> bH + followed by the decay H + —> t + v and its charge conjugate H~ process, 
and the search for high mass resonances decaying into tau-lepton pairs, the processes pp —> 
H/A — >• t + t~ . Both channels have been considered by the ATLAS and CMS collaborations 
and we briefly summarize below the resulting constraints on the /iMSSM. 

The CMS search [15] was performed with the 19.7 fb -1 data collected at y / s = 8 
TeV with the r-leptons decaying fully hadronically. 95%CL upper bounds have been set 
on the product of branching ratios BR(f— )• W7 + )xBR(L/ + — >t + u) from 1.2% at M H ± =80 
GeV (about the exclusion limit obtained on M H ± from LEP2 searches [18,19]) to 0.16% at 
M h ± = 160 GeV (beyond which phase space effects start to be too penalizing). The search 
excludes the entire range M H ± < 140 GeV for all values of 1 < tan /3 < 60. For larger 
masses, the areas where tan/3«8 at M H ± = 140 GeV to tan/3« 5-15 at M H ± = 160 GeV, 
in which the mt/ tan/3 component of the H±tb coupling is suppressed while the mb tan /3 
component is not yet enhanced, remain viable at the 95%CL. 

The ATLAS search for H± states [14] was also performed with the full 19.5 fb -1 data 
recorded at yfs = 8 TeV; the same channel as above, i.e. t —> bH + —> br + v with the tau- 
lepton decaying hadronically, has been used. Similar 95%CL upper limits than CMS have 
been obtained on the product BR(f->6iL + ) xBR(IL + — >t + v). Compared to the previous 
limits, a small additional area of the [tan/3, M H ±\ plane, at tan/3«8 for M H ± «100 GeV 
and to tan /3 ~= 6-10 at M H ± ~ 90 GeV, remains unexcluded by the ATLAS analysis, as 
a result of the presence of a large tt, W/Z+jet backgrounds in this mass bin. 

These limits can be turned into bounds in the [Ma, tan/3] parameter space assuming 
the usual relation M^± = + Myy. This is what is illustrated by the dark blue area 

of Fig. 11 in which the constraints on the [M^. tan /3] plan are shown: we take the limits 
on BR(t— > bH + ) x BR(i7 + —> t + u), that we calculate using the program SDECAY [68], and 
interpret them in the /iMSSM. We see that in our case, the entire area Ma ;$ 140 GeV 
is excluded at low tan/3 values, reducing to Ma < 130 GeV at high tan/3. In contrast to 
ATLAS and CMS, which use different means to calculate the product of branching ratios, 
we do not have the holes at the extreme values of M H ±. 

In the previous references, both experiments performed also searches for heavier H ^ 
states with M H ± >mt+mb~ 180 GeV, by considering the process pp-$H±t{b) with again 
H + ^t + v. The areas tan/3>45 (60) at M h ±k, 200 (250) GeV are excluded by the CMS 
collaboration (the ATLAS group assumed BR (H + — > t + v) = 100% in this area while it 
should be only 10%). These limits are much less powerful than those obtained from the tt 
search as will be seen shortly. However, there was in interesting CMS search in the channel 
H + —> tb performed with 19.7 fb" 1 data at y/s = 8 TeV [56]. Limits on a(pp —> tbH ± ) 
assuming BR(L7 + — > tb) = 1 have been set and surprisingly, one is not far from being 
sensitive to the very low tan (5 area and with about a factor of two more data, one would 
have probed tan/3 « 1 for M H ± « 200 GeV. 
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Figure 11: The [tan/3,M^] plane in which thepp— >H/A — > t + t~ (light blue) and t—>bH + —> brv 
(dark blue) search limits of ATLAS and CMS are combined and interpreted within the 6MSSM. 
The red-solid (blue-dotted) line is for the ATLAS (CMS) observed limit only in the M™ ax scenario 
and the region at the left of the black-solid line is the one where the 6MSSM is ill-defined. 

We note that these limits from searches exclude a substantial part of the low 
[tan ( 3 , Ma\ area in which the /iMSSM is ill defined, i.e. the region at the left of the thick 
black solid line in Fig. 11. The exclusion is valid in this area as all the ingredients used 
to obtain the limits depend only on M H ± and tan/3 and do not involve the CP-even 
Higgs parameters Mh or a that were undefined. The only assumption is that the relation 
Mjj ± = + Myy remains valid. Hence, one problematic issue within the 6 MSSM is in a 

sense partly solved by these “model-independent” exclusion limits (we do not address here 
the theoretical issue of the validity of the entire model at tan /3 values close to unity). 

The most important search mode in the MSSM is certainly the pp —>• H/A —>• rr 
channel. The ATLAS collaboration has searched for this signal using the 19.5-20.3 fb ” 1 
data collected at 8 TeV [16] while CMS has used the full 24.6 fb -1 data collected at 7+8 TeV 
[17]. Both collaborations consider the leptonic (+ 7 / 1 ), semi-leptonic (r^Thad) and hadronic 
(Thadihad) t decays and CMS also considers the case where an additional 6 -quark is present 
in the final state. Limits at the 95%CL on a(pp —> t + t~ ) as a function of the invariant 
mass M tt of the tau-lepton pair have been given by the two experiments. 

Our procedure to interpret these limits in the 6 .MSSM is as follows. First, we com¬ 
bine the ATLAS and CMS 95% CL exclusion limits on a(pp —> t + t~). We then compare 
them with the numbers that we obtain for the rates, namely the H/A production cross 
sections calculated with the programs HIGLU [52] and SusHi [51] and the branching ratios 
BR (H/A — > tt ) calculated using HDECAY [41]. These are derived assuming that the pa¬ 
rameters Mh and a given by eqs. (2.6) and ignoring all possible SUSY effects. For a given 
[tan/3, Myi] point, if in a window of A M rr = 40 GeV centered on M 4 we have either the 
observed h boson or the heavier H state, we include their contribution to the signal. 
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The result is shown Fig. 11 by the light blue areas of the [M^tan/?] plane and is 
truly impressive. The largest area is the one that excludes all values above tan /3 ~ 8 for 
Ma ^ 300 GeV, extending to tan f3 « 20 (40) for Ma < 700 (900) GeV. 

The ATLAS and CMS observed limits, when interpreted in the M™ ax benchmark 
scenario, are also displayed in Fig. 11. We observe that below Ma ~ 200 GeV, our limit 
is less restrictive, the reason being that we do not make use of any refinements in order 
to treat the regions in which the three Higgs bosons have a comparable mass and, also, to 
deal with the observed signal of the lighter h boson. In contrast, our limit is stronger at 
masses above Ma ~ 300 GeV, the reason being that while in our case there is no SUSY 
decays of the Higgs bosons, the M™ ax scenario leads to H/A decays into charginos and 
neutralinos that are significant. Indeed, for the gaugino and higgsino mass parameters of 
the M™ ax benchmark (and the slightly modified ones), /i = M 2 = 2M\ = 200 GeV, the LSP 
has a mass of m v o ~ 100 GeV while the heavier neutralinos and the charginos have masses 

A-l 

of the order of ^200 GeV, so that many SUSY decays occur starting from Ma~ 300 GeV 
and all of them will be present for Ma ~ 500 GeV and above. These decays will have 
substantial branching fractions, in particular in the intermediate tan/3 ~ 10 range where 
they become dominant, as can be seen from Fig. 9. Hence, the H/A —> tt branching ratios 
are suppressed in this this case, resulting in a weaker exclusion limit (which is unfortunate 
since this parameter configuration is almost certainly excluded by the direct searches [62]). 

A second step is to extrapolate the ATLAS and CMS limits to the low tan/3 region, 
which as discussed earlier, can be described within the hMSSM approach in contrast to 
the benchmark scenarios used by the collaborations. Two islands were discovered during 
the exploration. A first and substantial area is at very low tan/3 and Ma, tan/3 < 2 
and Ma < 230 GeV. Here, because part of the area is ill defined, we consider only the 
production and the decay of the A state that depend only on tan/3 and Ma as the decay 
channel A —> hZ , which introduces a dependence on the angle a through the AhZ coupling, 
is not yet kinematically open. In this domain, as discussed earlier, both the gg —> A 
cross section (dominated by the top-quark loop) and the branching ratio for the decay 
A — > tt (which, together with the one for bb and cc, is the only significant one to occur) is 
substantial. Hence, despite of the fact that we have only one resonance, the cross section 
times branching fraction is large enough to generate an observable signal. The excluded 
area from this search removes the small residual part of the ill-defined /iMSSM region, that 
was left after imposing the exclusion limit from the H± searches discussed above. 

More surprising at first sight, we discovered another smaller island at Ma ~ 350 GeV 
and tan /3 ~ 2-4. It turns out that, around the tt threshold, there is a very strong increase 
of the gg -A A amplitude as the form-factor Re(A]y 2 ) ~ 5 is maximal at the 2 mt threshold 7 , 
r = M\/Am/ ~ 1 . At the same time, BR(A — > tt) ~ is substantial being a few percent, as 
the other decays except for bb are slightly suppressed, A —>• hZ by the coupling ss cos(/3 — a) 
and A —> tt by phase space effects (only the three-body decay channel is kinematically open 
and it is suppressed). Hence, there would have been a substantial surplus of events from 
the gg—>-A—)■ tt process in this limited area that is excluded by the search. 

'In fact the NLO QCD corrections to a(gg^>-A) introduce a Coulomb singularity exactly at threshold 
[46]. However, we have checked that the rate increase in this new area is not due to this unphysical feature. 
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4.2 Interpretation of the bosonic Higgs decay modes in the hMSSM 

We now turn to the constraints that can be imposed on the [tan /3, Ma] plane by considering 
the bosonic decay channels of the heavier H and A states. In contrast to the H/A—^tt and 
H± —> tv searches, no interpretation of these modes has been done in the context of the 
MSSM by the ATLAS and CMS collaborations. In the following, we will therefore adapt 
the constraints that have been obtained either in the context of the SM but with a heavier 
Higgs state than the observed one, or in extensions of the SM other than the MSSM. We 
will focus on the experimental analyses that provide the most stringent constraints. 

The massive gauge boson channels H —> WW, ZZ have not been discussed in the 
context of the MSSM but important information can be borrowed from those performed for 
a heavy SM Higgs boson. For instance, a search of the H — >WW —±f,tvv and H —> ZZ— 
channels have been conducted in the SM with the full event sample recorded at the LHC 
first phase, i.e. 4.9 fb -1 at y/s = 7 TeV and 19.4 fb ” 1 at y/s = 8 TeV for H — > WW* 
[23] and 5.1 fb -1 at y/s = 7 TeV and 19.7 fb ” 1 at y/s = 8 TeV for H —> ZZ* [24], The 
high mass range was analyzed and the events corresponding to the observed state with a 
mass of 125 GeV were considered as a background. In the ZZ channel, an additional CMS 
analysis in the H —> ZZ* —> 2£2q channel has been made with the 19.6 fb -1 data collected 
at 8 TeV [69]. All these analyses exclude a significant area in the hMSSM parameter space 
at low and moderate tan /3. 

In the case of the H —> WW —> Mvv CMS search, when all production channels are 
included (there is a dominance of the gg —> H mode of course) and the various final state 
topologies are summed up, a Higgs particle with a SM-like coupling to gauge bosons is 
excluded from Mu ~ 200 GeV to « 600 GeV. The 95%CL upper limit as a function of 
Mu and relative to the SM expectation can be easily turned into an exclusion area in the 
[Ma, tan/3] plane by considering the production and decay rates of the MSSM H state 
discussed in section 3. The result is shown in Fig. 12 where the area excluded by this 
search, interpreted in the context of the hMSSM, is depicted in dark green. 

The exclusion area starts at relatively high tan/3 values, tan/3> 10, and light A, 
Ma« 140 GeV (below this limit, we enter the domain in which the model is ill defined, a 
domain that extends to Ma ~ 160 GeV and tan/3 ~ 1) where one has an H state with a 
mass Mu > 160 GeV and a coupling guw = cos(/3 — a) that is not very small as shown in 
the right-hand side of Fig. 3, allowing for substantial H production times decay rates. For 
tan/3~l, the excluded region extends to M,4~250 GeV, when other decay channels such 
as H^-hh and even H —ttt open up and suppress the massive gauge boson decay modes. 

As a result of its clean final state and despite of the low statistics, the H—>ZZ— 
search turns out to be more constraining at high mass and excludes a SM-like Higgs boson 
up to Mu ~ 800 GeV (with a search domain extending to 1 TeV). While for low Mu 
values, H —> ZZ is less powerful than the companion H — >WW mode as a consequence of 
the reduced phase space, it clearly becomes the leading channel for Mu </250 GeV. In fact, 
because of the higher statistics, the most severe constraint is obtained in the combination of 
the H—> ZZ—>4£, 2£2v, 2£2q topologies that was performed in Ref. [69]. Here, the 95%CL 
exclusion of a Higgs state with SM-like couplings extends to a mass close to 1 TeV. 


- 30 - 


The area excluded at 95%CL by the non observation of these ZZ final states at the 
LHC outside the M/, ~ 125 GeV mass window is given by the light green area of Fig. 12. 
It extends from Ma ~ 160 GeV to Mh~ 280 GeV and concerns all values tan/3 < 5. An 
additional small area around Ma ~ 300 GeV and tan /? < 2 , in which Mh is close to the 
2 mt threshold and the gluon-fusion amplitude is maximal thus enhancing the gg^H 

cross section, is also excluded. 

One should note that in a dedicated MSSM search, not only this H —> ZZ channel 
but also the H —> WW mode will lead to more effective constraints as the SM and MSSM 
Higgs particles have total decay widths that are completely different at high masses [5]. 
Indeed, while the SM state would have been a very wide resonance, the MSSM H boson 
is a relatively narrow resonance as shown in Fig. 5, allowing to select smaller bins for the 
VV invariant masses that lead to a more effective suppression of the various backgrounds. 



Figure 12: Constraints in the [tan/3,ALi] plane of the 6 MSSM from search at the LHC for a 
heavy CP-even boson decaying into WW, ZZ, hh and a heavy CP-odd boson decaying into hZ 
final states. The searches are for ^fs = 7 + 8 TeV c.m. energy and 25 fb -1 of accumulated data. 
The dashed area is the one that is ill-defined in the 6 MSSM. 

The resonant H -^hh channel, which is important in the mass range between 250 and 
slightly above 350 GeV has been considered by both the ATLAS and CMS collaborations 
with the ~ 20 fb -1 of data collected at -y/s = 8 TeV. The main focus was on the 77 56 
signature [25, 26] but additional searches in the 4 6 -quark final state have been recently 
reported [70,71]. However, neither collaborations has interpreted the 95%CL exclusion 
limits in these channels in the context of the MSSM, the main reason being again that 
the low tan (3 area in which these signals occur is not theoretically accessible in the usual 
benchmark scenarios used for the MSSM Higgs sector. 

The interpretation is however straightforward in the LMSSM as the trilinear self¬ 
coupling A Hhh that controls the H -+ hh decay rate is simply given, as shown in eq. (2.10), 
in terms of the angles a, /3 and the radiative correction matrix element AM | 2 that is fixed 
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in terms of tan/3 and Ma if the constraint Mh = 125 GeV is used. We have adapted the 
constraints from these analyses to the LMSSM case and the resulting excluded domain in 
the [tan 0,Ma\ plane is shown in purple in Fig. 12. It covers the very low tan/3 region, 
tan/3 < 2 , for the mass range between Ma ~ 270 GeV (which implies Mu > 250 GeV for 
these low tan/3 values) and Ma ~ 330 GeV, i.e. slightly before the 2 mt threshold. 

Similarly to the previous channel, the A —> hZ mode has only been considered in the 
context of two Higgs doublet models [37] and not in the MSSM. A CMS analysis considered 
the final state bb£ + £~ with the ~ 20 fb _1 collected in 2012 at 8 TeV [27]. A search of both 
the A —>• hZ and H —>• hh channels has been performed by CMS again in the multi-lepton 
and eventually photon finale states [72]. The impact of the 95%CL exclusion limits of these 
studies, when interpreted in the context of the /iMSSM, is illustrated by the yellow area 
of Fig. 12. The ranges tan/3 < 3 and Ma ~ 230-350 GeV should be in principle excluded 
with the present data. 

4.3 Summary of the constraints at 8 TeV and projections for 14 TeV 

Wrapping up the discussion up to this point, the impact on the [Ma, tan/3] plane of the 
searches in the fermionic Higgs decays H/A —> rr and H± —> tv and in the bosonic ones 
H — > WW, ZZ, hh and A —>■ Zh performed by the ATLAS and CMS collaborations at 
t/s = 7+8 TeV with up to ~ 25 fb _1 data are combined in Fig. 13. The outcome is very 
impressive. The high tan /3 > 10 is entirely excluded for Ma 500 GeV by the rr searches. 
The range tan/3 < 4 is excluded for Ma < 250 GeV by the H -+ WW, ZZ channels. For 
tan/3 < 2, the excluded domain extends to Ma < 350 GeV when the channels H —> hh 
and A —> hZ are considered. The entire low Ma region, Ma < 140 GeV, is excluded by 
the H + —y tv search for any value of tan /3. An additional portion of this low Ma area is 
excluded by the A -+ rr search for low tan 0 values. 



Figure 13: The combined constraints in the [tan /?, Ma] plane of the hMSSM from searches at the 
LHC for the heavier H, A and H ± bosons decaying into either fermionic or bosonic final states. All 
the searches performed at c.m. energies up to ,/s = 8 TeV and 25 fb -1 data are included. 
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In fact, the entire area in which the /rMSSM is not mathematically defined, and which 
is delineated by the solid line in the figure, is excluded by these H± and A searches that 
do not involve the undefined CP-even H boson mass Mh and the mixing angle a. 

These constraints, if no new signal is observed, can be vastly improved at the next 
phase of the LHC with a center of mass energy up to y's = 14 TeV and with one or 
two orders of magnitude accumulated data. More optimistically, this implies that the 2er 
sensitivity for a heavier MSSM Higgs boson will be drastically enhanced at the next LHC 
phase. Starting from the expected median 95%CL exclusion limits that have been given by 
the ATLAS and CMS collaborations in the various searches performed at 8 TeV with ~ 20 
fb -1 , we have made an extrapolation to this next LHC phase with yfs = 14 TeV and 300 
fb -1 data. We have naively assumed that the sensitivity will simply scale with the square 
root of the number of expected events and did not include any additional systematical 
effect. This comes from the observation that the results of the experimental analyses are 
limits on the signal cross section at a given c.m. energy for a given resonance mass bin, 
for a channel that is subject to a given background rate R^(Ma) at this mass 
bin, when the integrated luminosity is fixed at a value £Knowing the sensitivity limit 
R§ ( Ma ) at = 8 TeV, one derives the associated limit at y / s = 14 TeV using 

RUm A ) = ^£ 8 /£i4 X ^Rf 4 {M A )/Ri{M A ) X Rl(M A ) (4.1) 

Having the knowledge of only the signal cross sections a s r (MA) for the various points and 
not the corresponding background rates, we assume that the latter simply and very naively 
scale like the signal cross sections. This is the case of some channels of interest, such as 
gg —> H / A —>■ tt whose main background is gg —> tt and as both are gg initiated processes, 
they roughly scale with the gg luminosity at higher energies. However, for many other 
channels such as H —> WW, ZZ or A/H —> tt, the irreducible background is mostly due 
to qq annihilation which increases more slowly with energy than the initiated gg signal 
processes. This makes our approach rather conservative. 

With this assumption, one obtains for the sensitivity at y / s = 14 TeV, R,/ a (Ma), needed 
to set the exclusion limit, that we turn into a 95%CL sensitivity, for a given Ma 

Ru(Ma) « VAACu x ^af 4 (M A )/al(M A ) x Rl(M A ) (4.2) 

The output of this procedure is presented in the [tan/3, M A ] hMSSM plane in Figs. 14 
for the fermionic (left) and bosonic (right) Higgs search channels. In the former case, we 
have included in addition the channel pp —> tbH + —> tbtb which now shows some sensitivity 
a low tan/3 and not too high M H ± values. The combined expected 95%CL sensitivities 
are shown in Fig. 15 and, as can be seen, a vast improvement of the current sensitivity 
to the MSSM parameter space is foreseen in all channels. This is particularly the case of 
the A/H -A tt channels which alone, closes the entire region below M A 350 GeV for 
any tan /3 value, while the H -A WW, ZZ modes which show sensitivity up to M A ~ 600 
GeV at very low tan/3. In the Higgs mass range in which they are relevant, i.e. below the 
tt threshold, the channels H -A hh and A —> hZ start to probe rather high tan/3 values, 
tan/3 ~ 4 and tan/3 ~ 6 , respectively. 
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Figure 14: Projections for the LHC with y/s = 14 TeV and 300 fb -1 data for the 2cr sensitivity in 
the /iMSSM [tan/3, Ma] plane from the search for A/H ± states in their fermionic decays (left) and 
A/H states in their bosonic decays (right). The same color code as at yfs = 8 TeV has been used 
and, for the fermionic channels, we add a constraint from the H + —> tb mode depicted in dark blue. 



A/H —>TT 
H^-rv 
rf^tb 
H -> WW 

H->ZZ 

A-> Zh 

H -> hh 


Figure 15: Projections for the LHC with s = 14 TeV and 300 fb -1 data for the 2cr sensitivity 
in the /iMSSM [tan/3, Ma) plane when the searches for the A/H/H ± states in their fermionic and 
bosonic decays are combined. 


Nevertheless, there will remain an area of the /iMSSM parameter space, at tan/3 < 4 
and masses above Ma ~ 400 GeV to name it, which will not be accessible by the channels 
that have been considered so far in the search of the heavier H/A and states. To probe 

this area, the high luminosity option of the LHC with C = 3 ab _1 data or a higher energy 
pp collider, such as the presently discussed Fcc-pp at \fs ~ 100 TeV will be necessary. 
However, as it was discussed in many instances in this paper, this virgin area is the ideal 
territory to perform searches in the gg —>■ H/A —> tt channel to which we turn our attention 
now. 
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4.4 Including the pp —>■ H/A —> tt channel 

As it was discussed at length in the previous section, for low tan (3 and high Ma values, 
the decay modes H/'A tt of the heavier MSSM scalar and pseudoscalar Higgs states will 
largely become the dominant ones while the gg —> H/A cross sections are still substantial 
thanks to the large Higgs coupling to the top quarks that mediate the production process. 
Hence, the search for resonances decaying into tt final states will be mandatory in order to 
probe these areas of the [Ma, tan/3] parameter space at the LHC. However, a peak in the 
invariant mass distribution of the tt system, that one generally expects in the narrow-width 
approximation, is not the only signature of a Higgs resonance in this case. Indeed, the gg —> 
H/A signal will interfere with the QCD tt. background which, at LHC energies, is mainly 
generated by the gluon-fusion channel, gg —>• tt. The signal-background interference will 
depend on the CP nature of the d> = H/A boson and on its mass and total decay width; 
it can be either constructive or destructive, leading to a rather complex signature with a 
peak-dip structure of the tt invariant mass distribution. 

These aspects are known since already some time and have especially been discussed 
in the context of a heavy SM Higgs state [73] and, hence, for the CP-even Higgs case. 
The slightly more involved MSSM situation, as there are one CP-even and one CP-odd 
resonances that are close in mass, has been addressed only in a very few places; see for 
instance Refs. [74,75]. Dedicated analyses have been performed at the parton-level only 
and do not make use of recent developments like boosted heavy quark techniques [76] 
that could allow to enhance the observability of the Higgs signal. The ATLAS and CMS 
collaborations have performed searches for heavy states decaying into tt pairs [77, 78] but 
did not specifically address the complicated Higgs situation as only electroweak spin-one 
resonances, like new neutral gauge bosons or electroweak Kaluza-Klein excitations, were 
considered. In these two cases, the main production channel is qq annihilation and there 
is no interference with the (colored) QCD qq —> tt background and the resonance signal 
simply appears as a peak in the invariant mass distribution of the tt pair. 

A full and realistic Monte-Carlo simulation of the gg —* H/A —»• tt process including 
the effects of the interference and taking into account reconstruction and detector aspects 
is beyond the scope of this paper, and will be postponed to a future publication [79]. Here, 
we will simply make a very crude estimate of the sensitivity that can be achieved in this 
channel, relying on previous ATLAS [77] and CMS [78] analyses performed at yfs = 8 TeV 
c.nr. energy in the spin-one resonance context mentioned above. We will naively consider 
the number of signal and background events, applying very simple kinematical cuts and 
ignoring the complicated interference effects, and delineate the area in the [tan /3, Ma] 
/iMSSM parameter space in which one has AQ gna i/i^/Wkg > s. The significance s = 5 
would correspond to a 5cx observation of the Higgs signal while s = 2 would be a first 
hint of the new effect; in the absence of any effect, s = 2 would correspond to the 95%CL 
exclusion limit of the phenomenon. To further simplify our analysis, we will assume that 
the two heavy A and H states are mass degenerate so that the signal rate is simply the 
sum of the A and H production cross section times the respective branching ratios in their 
decays into tt pairs (which, as we have already seen, is a good approximation). 
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The main ingredients of the analysis are as follows. The normalization of the Higgs 
signal has been obtained using the programs HIGLU for the production cross sections and 
HDECAY for the decay branching ratios. The total cross section of the SM background 
(which will serve as a normalization) has been obtained using the program Top++ [80]. For 
the input rrit = 173.2 GeV one obtains for the background rate at the first stage of the 
LHC with y/s = 8 TeV 

o^ D (PP ~>tt) = 247.7 tgl +^5 pb at yfs = 8 TeV (4.3) 

when the renormalisation and factorization scales are fixed to pr = Hf = n r h- In this 
equation, the first error is the one due to the scale variation within a factor of two from 
the central scale, and the second one the PDF+a s uncertainty. This value for the cross 
section is obtained at NNLO in QCD including the resummation of next-to-next-to-leading 
logarithmic (NNLL) soft gluon terms and it turns out that it is only 3% larger than the 
value of the cross section when evaluated at NNLO [82], Note that at yfs = 14 TeV, using 
the same approximation and ingredients, one would obtain for the cross section 

*2? D (pp tt) = 966.0 till - 4 oJ Pb at yfs = 14 TeV (4.4) 

Using the program MadGraph5 [83], we have generated the signal and background 
cross sections for the process pp —> tt. The differential cross section as a function of the 
invariant mass of the tt system, da/dmtt, is shown at yfs = 8 TeV in the upper part of 
Fig. 16 where mass bins of 10 GeV have been assumed. We overlay on the continuum QCD 
background distribution (in black solid line), the distributions for the A signal only (the 
colored lines) with tan f3 = 1 and three possible mass values, Ma = 400, 600 and 800 GeV. 
In order to see the signals in the figure, we have multiplied the distributions by a factor of 
5, 50 and 300, respectively. 

In order to enhance the significance s, one could apply very basics kinematical cuts 
that suppress the background while leaving the signal almost unaffected. In the left and 
right-hand sides of Fig. 16, we show two distributions (as we are interested in the shapes 
only, the distributions have not been re-weighted with the correct A"-factors etc., and 
the integrated areas thus correspond to the Monte-Carlo cross sections). The first one is 
1/a x da/d cos 9* where 9* the helicity angle between the off-shell Higgs boson boosted 
back into the top quark pair rest frame and the top quark pair direction (left). As can be 
seen, while the signal distribution is almost flat, the background is peaked in the forward 
and backward directions; a cut | cos 9 * | < 0.8 for instance would remove a large sample 
of background events. A second distribution per 10 GeV bin is in terms of the transverse 
momentum of the top quarks, 1 /ax da/dpx (right). They show a characteristic behavior 
for the signal events, with a pronounced peak and then a sharp drop. One grossly estimates 
that, for the mass value Ma = 800 GeV for instance, a cut on the p l T distribution could 
allow to suppress the background by a factor of « 6 . 

Assuming that when applying all kinematical cuts, one could suppress the tt QCD 
background by an order of magnitude without significantly altering the Higgs signal, 
we delineate in Fig. 17 the regions of the [tan /3, Ma] plane in which one would expect 
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Figure 16: The pp —> tt signal and background at y/s = 14 TeV for Ma = 400, 600 and 800 GeV 
and tan j3 = 1: the invariant nitt (top), cos 9* (bottom left) and p l T (bottom right) distributions. 

-^signal/y^Vbkg > 2,3,4, 5. For the previous LHC run with y/s = 8 TeV and 25 fb ” 1 data 
(left), one observes that a two 2 a “evidence” for a new resonance, or a 95%CL exclusion 
of the relevant point of the MSSM parameter space in the absence of any effect, can be 
achieved for tan /3 ~ 2.5 if Ma ~ 350 GeV and tan /3 ~ 1 when Ma ~ 550 GeV. A 5er 
discovery could be achieved in this case in a much smaller area of the parameter space, 
tan/3 ~ 1.5-1 and Ma ~ 350-450 GeV. The situation could be vastly improved at the next 
LHC run with yfs = 14 TeV and 300 fb ” 1 data (right) as one could be sensitive to tan/3 
values tan/3 « 7 for Ma ~ 350 GeV and mass values Ma ~ 1 TeV for tan/3 « 1. 




Figure 17: Sensitivity in the gg — >• H/A —>• tt channel in the [tan/9, Ma] plane of the hMSSM at 
the 2,3,4, 5cr level y/s = 8 TeV and 25 fb -1 (left) and yfs = 14 TeV and 300 fb -1 (right). 

The 2cr sensitivity in the [tan /3, Ma] plane when the H/A —> tt reach (using the the 
assumptions above) is superimposed to the sensitivity in all the fermionic and bosonic 
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channels discussed previously is displayed in Figs. 18 and 19 for, respectively, the previous 
and the next LHC phases. As can be seen, a vast improvement in the sensitivity is expected 
if the H/A—ttt channel is included, in particular at the forthcoming LHC run with y/s = 14 
TeV and 300 fb ” 1 data. The improvement is even more impressive at the high-luminosity 
LHC option, when the luminosity is increased to 3000 fb -1 ; see Fig. 20. In this case, almost 
the entire /iMSSM parameter space, up to Ma values close to ~ ITeV, can be probed. 



Figure 18: Expectations for the 2a sensitivity in the /iMSSM [tan 8, Ma] plane when the searches 
for the A/H/H ± states in all channels, including the gg —> H/A —► tt process, are combined at the 
LHC with y/s = 8 TeV and 25 fb -1 data. 



Figure 19: The same as in Fig. 19 but at the LHC with </s = 14 TeV and 300 fb 1 data. 
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Figure 20: The same as Fig. 20 but at the high -luminosity LHC option with 3000 fb 1 data. 

5. Conclusions 

In this paper, we have addressed the issue of covering the entire parameter space of the 
MSSM Higgs sector at the LHC by considering the search of the heavier H , A and 
states that are predicted in the model, in addition to the already observed lightest h boson. 
These searches should not only be restricted to the channels that have been considered so 
far by the ATLAS and CMS collaborations, namely those with a surplus of tv events 
and those with high mass resonances decaying into r lepton pairs, which would signal the 
presence of new contributions from the t —> bH + —> brv and pp —>• H/A rr processes, 

which are mainly relevant for the high tan /3 region of the MSSM Higgs sector. Search for 
heavier Higgs bosons should also be conducted in channels that are more appropriate for 
the probing of the low tan (3 region and which, until now, have been overlooked. 

We have first discussed and refined the /iMSSM approach introduced in Ref. [10] in 
which the dominant radiative corrections to the MSSM Higgs sector, that introduce a 
dependence on numerous SUSY parameters, are traded against the measured mass Mh = 
125 GeV of the Higgs boson which was observed at the LHC, thus allowing to describe 
again the entire Higgs sector of the model with only two input parameters. This simple, 
economical and “model independent” approach permits to reopen the low tan (3 region, at 
the expense of considering the possibility that the scale of SUSY-breaking is extremely 
high, Mg 1 TeV, and that the model is severely fine-tuned. The /iMSSM is expected to 
be viable down to values tan/3 ~ 2 and, for higher tan/3 values, reproduces to a very good 
approximation the standard results of the MSSM Higgs sector. This is particularly true 
if the higgsino mass parameter is much smaller than the SUSY-breaking scale, p, <C Ms, 
an assumption that is natural at low tan (3 values which imply a very high SUSY-breaking 
scale. Thus, searches for new signals in the MSSM Higgs sector can be performed in the 
entire [tan/3, M^] parameter space, in a reliable way for tan/3 > 2. Nevertheless, in an 
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effective approach, one can eventually extrapolate to values of tan/3 very close to unity, 
despite of the fact that the scale Ms required to reach this value is so high that its renders 
the model not only too fine-tuned but also potentially inconsistent. 

We have then analyzed the production and decay modes of the H , A and H ± particles 
at the LHC, with a special attention to the low tan /3 region in which the top quark plays a 
prominent role, as its couplings to the Higgs bosons are not strongly suppressed compared 
to the SM case. We have first shown that the searches that are presently conducted by 
ATLAS and CMS can also be relevant at low tan /3. This is for instance the case of the 
pp —> A —> tt and pp — > tbH + — > tbtb processes at low to moderate Ma values. We have 
then shown that search channels such as H WW, ZZ, hh and A —>• hZ, when interpreted 
in the context of the ZiMSSM, can also probe the low tan /3 and not too high Ma regions. 
In fact, already with the 25 fb -1 data collected at energies up to ^fs = 8 TeV, the ATLAS 
and CMS collaborations exclude the possibility tan /3 < 4 ( 2 ) up to pseudoscalar Higgs 
masses of Ma ~ 250 (350) GeV. At the upcoming stage of the LHC, with an expected 
energy and luminosity of s/s = 14 TeV and £ = 300 fb _1 , the entire parameter space i.e. 
for any value of tan j3 could be probed up to Ma ~ 400 GeV, when combining the searches 
in the usual fermionic channels and in the additional bosonic channels discussed here. 

An important message conveyed by the present paper is that, in order to fill or close the 
gap in the MSSM [tan f3, Ma] plane left by the fermionic and bosonic searches mentioned 
above, one should definitely consider the pp —>• H/A —>• tt process. Indeed, at low tan/3 
and for Higgs masses above the tt kinematical threshold, the decays H/A —> tt become 
the dominant ones, suppressing the rates for the other decay channels to a very low if not 
negligible level. On the other hand, the gg —> H/A production mode has a still significant 
cross section as the top quark that generates this loop process has substantial couplings 
to the H/A states at sufficiently low tan/3 values. This is not a very easy search channel 
in view of the formidable pp —> tt QCD background. Nevertheless, it exhibits very special 
and interesting features such as an interference with the QCD background that leads to a 
rather involved peak-dip structure of the signal. 

We have not performed a detailed and realistic study of this process but attempted to 
roughly quantify the observation of a signal at the LHC, relying on present ATLAS and 
CMS analyses in searches for heavy (non Higgs) resonances decaying into top quark pairs 
at 8 TeV center of mass energies, and discussed its possible implications. It appears that 
the channel gg —>• H/A —>• tt, would be capable of covering partly the area at low tan/3 
and high Ma, hence allowing for a full coverage of the [tan /3,Ma\ plane of the MSSM 
up to Higgs masses Ma ~ 600 GeV with 300 fb -1 data at yC = 14 TeV. At the high 
luminosity option of the LHC with 3000 fb -1 data, one could reach a full coverage of the 
MSSM parameter space for pseudoscalar masses closer to Ma ~ 1 TeV. 

More refined analyses are required in order to firmly establish the viability of the 
various processes discussed here, in particular the H/A —> tt channel. In view of the 
important role that it could play in the probing of the MSSM parameter space, the latter 
process is worth investigating in a more realistic way, including the interference between 
the Higgs signal and the QCD background. This is what we plan to do in a forthcoming 
publication [79]. 
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